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ABSTRACT 


The knowledge of wave refraction is important in many studies. The 
rapid and relatively easy gaining of this knowledge is made possible by 
the use of the modern high-speed digital computer. Large numbers of 
Spectral periods and incoming directions are easily investigated, and 
immediate results are obtained by use of the plot of the wave crest re- 
fraction from the computer. This program presents the wave crest refrac- 
tion pattern of numerous wave ray points rather than the single ray 
following technique. Its use is valuable in amphibious operation plan- 
ning, and in other studies where the distribution of wave energy along 


the shore is desired for the many periods of the wave spectrum. 
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la-Intseduetion 

The knowledge of the refraction pattern as an ocean wave approaches 
the shore is necessary in the determination of the energies of the wave 
in the surf zone. Wave energy is related in turn to many near shore 
processes or operations, such as beach erosion, sediment movement, or 
amphibious landing operations. 

While this investigation is slanted toward the needs for amphibious 
operations, the knowledge and methods used are amenable to any users who 
desire the pattern of wave crest or refraction, or an estimate of the 
ratio of the wave height at a point to the deep-water wave height (H/H,) 
for other purposes. 

Early work on the subject was concerned mostly with hand calculation 
of the various determining parameters and graphical method for construc- 
ting the refraction diagrams dao. Later, in early 1962, attempts were 
made to evaluate wave refraction using high-speed digital computers [ 3]. 
The most usable results published were those obtained by Harrison and 
Wilson [z]. While the work of Griswold and Nagle as well as that of Har- 
rison and Wilson gives worthwhile results, both studies are very limited 
with respect to the ratio of work input to the results obtained. They 
are concerned with following one ray at a time from deep water to the 
shoreline. To get a complete refraction diagram many inputs to the com- 
puter are required, and hand plotting of the results are necessary. 

The procedure outlined in this paper follows a wave crest composed 
of a number of ray points to the shoreline with immediately available 
results. No hand plotting or further program input is required for the 
rays with one period and the same initial deep-water direction. By 


changing only the wave period or wave direction, the spectrum of periods 


and range of angles can be investigated rapidly with little additional 
time involved. 

All of the computer programs for the investigation of wave refraction 
utilize the same basic procedure: following the wave ray (orthogonal) to 
the shoreline. The new program uses a depth field as an interpolating 
surface rather than the velocity field of (3) ձոզ [4]. Several interpo- 
lation surfaces have been investigated to represent the velocity or depth 
Field used (41. This program utilizes a quadric surface for interpola- 
tion of the depth field. 

The present study uses an input of deep-water wave period and wave 
direction to a computer program for determining wave refraction in a 
method such as proposed by Munk and Arthur, which employs the wave para- 
meters listed in Appendix tv [5]. An additional subroutine of the pro- 
gram computes the coefficient of refraction, K, and the ratio of wave 
heights, НАН: The latter two values are recorded along with the X and 
Y coordinates of each point along the wave crest. As a rapid means of 
viewing the refraction pattern of the wave crest, a graphical output is 
included which contours every third wave crest computed. Other parameters 


may be included in the output at the user's discretion. 


2. Method. 

The first step in the utilization of the program is the construction 
of a grid of depth values for the desired area. This grid must include 
starting points in deep water for all rays to be followed, such that the 
ratio of the depth to the deep-water wave length, d/Lo; is greater than 
0.5. Since the program is arranged to follow the waves from deep water 
to the shore, the grid origin must be in deep water. The convention is 
that the X-axis will be positive and increasing toward the shore while 
the positive Y-axis is 90 degrees to the left of the X-axis as shown in 
Figure 1. The grid interval is selected such that, in a given cell, the 
bottom contours are reasonably parallel to one another. 

Actual or interpolated depths at the grid intersections are recorded 
to the best accuracy available from the chart, and all actual depths are 
made positive. Extrapolated depth values are continued on land for two 
grid units from the shoreline and are made negative. For depths on the 
shoreline itself, zero is used. Depths on land outside of two grid units 
from the shoreline are made some arbitrary, negative, non-zero value. 

The procedure for assigning negative values for the nearshore land depths 
is required in the fitting of a surface to the localized depth values. 
The zero land depths are used for contouring the shoreline on the output 
graph. 

3. Example of Input. 

As an illustration of the procedure and to test the results, the 
southern portion of Monterey Bay, California, was used. The depth grid 
was selected so that the origin of the 18 rays would be positioned in 
water depths greater than 1,024 feet where але = 0.5 for a period of 20 


seconds. The direction toward which the rays proceed is 125° TRUE, making 


an angle with the X-axis (oriented at 90° TRUE) of - 35°. The grid inter- 
val selected was 1,500 feet, and the depth input was in fathoms to facil- 
itate the chart reading. The depth values were determined from the U.S. 
Coast and Geodetic Survey chart 5403. 

Other parameters that were required for input were the X and Y values 
of the starting point for the first ray. The values of the ray separa- 
tion parameter, B > are input as Bl and B2, and are always equal to 1.0 
in deep water. The period of the waves and the angle with respect to the 
X-axis are entered a$ T and Al, respectively. The following parameters 
are: NOR, the number of rays that are being followed; DIST, the distance 
between these rays; TIME, the time interval for advancement of the wave 
front; GRID, the grid interval; MM, the number of grid points in the X 
direction; and NN, the number of grid points in the Y direction. The 
angle of the wave direction is in relation to the X-axis and is the di- 
rection toward which the waves are moving. It can be positive or negative 
with respect to the X-axis. 

4. Computer Operations. 

The computer first reads the parameter values, then the depth grid 
as a column of Y-values for a constant X-value. The depths are converted 
to feet immediately. The subroutine DEPTFUN is called to compute the 
depth at the first point by fitting the closest nine grid point depths 
around the starting value to a quadric surface by the least squares meth- 


od, using an equation of the form: 
DEPTH = А, + А„Х + AsY + A XY + ДЕ” + Ag y^ 


where DEPTH is the value of the depth at that point, А, 2. аге constant 
9 e ~ 


coefficients, and X and Y are the distance values for the point. Each 


time a new depth is encountered the surface of best fit is computed from 
the surrounding nine grid points. Also included in the DEPTFUN subrou- 
tine is the computation of the wave velocity at that point for the given 
depth. This is done by an iteration process using the common result from 


solving the wave equation as used in H. 40; Pub. 234 [2]: 


Ը» ы tanh (204) ғ 


where C = wave velocity, g = acceleration due to gravity, T = wave per- 
iod, and d = water depth. As with other wave refraction investigations 
it is assumed that the wave velocity is a function only of water depth 
and wave period. Other factors such as bottom friction, percolation, 
currents, reflection, and wind are considered as not affecting the re- 
fracting waves. 

The wave ray is moved to the next point by solving for the ray curva- 
ture and projecting the ray forward in the time interval specified at the 
speed calculated for the point. 


The curvature is calculated using the expression [5]: 


ЄК = Z [sma f - cos4 gj 


where FK = ray curvature, and A = approach angle. 
To determine the values of X +1° Yael? А 817 апа PK 41 for the suc- 


ceeding point, an iteration process is used to solve the equations [3]: 
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Уы = Ул tD sin A 


where D = the incremented distance (D = CT) between points n and n+l, 
At the point n+l, the value of Beta is calculated by solving the 


second-order, non-linear differential equation [5]: 


DA + pee da 26 = О 


Ds? 
where 2 
ο - ւմ 295 
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The above equation is solved by the finite difference method 3 |. This 
results in an equation for the Beta value at the n+l point in terms of 
the Beta values at the two previous points. The equation to be solved 


is then: 





2 
/ а (Рр - 2) В, + (+- 245) β 
3 
2 + рр 
where p, q and D are defined above, and Ох and Ё, are the Beta values 
of the two previous points. It is to be noted that the calculation of 
Beta is made for the n+l point. In the BETA subroutine the coefficient 


of refraction is calculated by the relation: 


K = 1/6 


where K = the refraction coefficient. 
To determine the ratio of wave height at the point to the deep- 


water wave height, the shoaling factor has been approximated for values 


10 


of с/с, < ] by a curve such that [2] : 


= В £ с 2 eve 
Hs zc ro /2.8150 (7) + 28.88% (£) τ 29 9257 (2) 


4 
СА _//. 28/5(6 ) 


As Co» the deep-water wave velocity, is a function of period alone 
and C is known for any point (X,Y), the equation may be evaluated at any 
or every point at which the calculations are made for C. The ratio of 


wave heights ԳՆԴԸ) is found from the equation [2]: 


H/H = Hs 1/8 


5. Output from Computer. 

The printed output from the computer consists of the X and Y values 
of each position after advancement by the increment of time specified. 
This is in units of yards for ease in hand plotting on the charts. The 
coefficient of refraction and the wave height ratio for each point are 
given. 

Included in the program output is a graphical plot of the wave crests, 
which is programmed for the Calcomp 160 computer system, utilizing the 
DRAW subroutine of Appendix II. The first and every third crest there- 
after are plotted and points connected depicting the wave crest. The 
scale is such that a true representation of the wave crest is presented. 
From this graph and if desired, hand plotting, areas of convergence and 
divergence are easily seen. By knowing the number of the crest, the para- 
meters of refraction coefficient and wave height ratio can be found from 
the printed data. 

6. Program Development. 


To calculate the coefficient of refraction and the wave height ratio 


ae 


it is necessary that the interpolation surface for the depth values be of 
the second order, at least, so that the second partial derivatives would 
be available for the computation of these values. Consequently, a quadric 
surface was used for the representation of the depth values at the closest 
nine grid points to the position at which the values were desired. The 
partial derivatives of the surface in the X and Y directions were used to 
evaluate the partial derivatives of the velocity function as proposed by 
Harrison and Wilson [4]. However, a power series representation of the 
hyperbolic tangent was used to evaluate the velocity derivatives rather 
than the method used by Harrison and Wilson, as shown in Appendix III. 

The calculation of the Beta parameter is required for determining 
the refraction coefficient. The finite difference method is used in this 
program. However, the use of this method requires that the distances be- 
tween the points at which the equation is being solved be equal. In the 
method presented here where the distance is a function of the velocity at 
each point, this does not hold exactly true. The difference in the veloci- 
ties at the successive points is on the order of one foot per second due 
to the shallow contour gradient in the particular area of interest. This, 
of course, will not hold true for all cases. A better method would be to 
solve the second-order equation in Beta by the Kelvin method which re- 
quires only the distance between the point n and n+l, which is readily 
obtainable. 

7. Discussion of Results. 

The printed results are shown in Appendix I for several of the wave 
crests. The X and Y values from this type of presentation were plotted 
by hand as in Figure 1. This plot shows areas of marked divergence along 


the shore. When this figure is compared with Figure 2, the result of the 


12 


graphical method of H.O. Pub. 234, little difference is noted suggesting 
that both the computer method and the graphical method produce similar 
results. This is the aim of the investigation. Other comparisons were 
made using different directions and periods with comparable results. 

In Figure 1, the seventh and eighth wave rays are seen to cross and 
continue to the shore. This crossing is attributed to the bottom contours 
of ray eight when the ray first reaches shallow water (d/Lo 20.5). There 
is a small but steep-sided indentation in the otherwise gradual slope of 
the area. This causes the ray to bend due to the steep depth gradient 
which results in convergence with ray seven below it. From this point 
to the shore the refraction is similar to that of ray six. Figure 3, the 
computer drawn plot, shows all of the wave crests from the first advance- 
ment to the shore. The ray crossing is evident. 

Ihe values of the coefficient of refraction and the wave height ratio, 
Appendix I, show little refraction for crests one and two as expected, 
since the waves for the most part are in deep water. However, later values 
of these parameters do show the effects of the refraction seen in Figures 
1 and 3. The values of these parameters require verification. The values 
do appear to be qualitatively reasonable for the refraction encountered, 
and when compared to the values estimated by the technique of measuring 
the distance between orthogonals. 

It is estimated that the time necessary to construct the grid and 
record the depth values is three to four hours depending on the size of 
the area desired and the gradient of the bottom contours in that area. 

A shallower gradient requires a larger area to ensure that the wave rays 
start in deep water. The time required to punch the data cards and to 


check the results is a function of the experience of the operator. The 


15 


computer time is approximately five minutes for the computation of the re- 
sults in this paper, Appendix I, plus an additional five minutes to draw 
the graph. This time can be reduced drastically if the program as com- 
piled by the computer is put on tape and input in that form. The compiler 
takes over half of the computer time at present. The time is spent con- 
structing the grid of depth values only once. From this grid all of the 
various wave periods and directions can be investigated without further 
effort of constructing another grid as is necessary in the hand plotting 
of the refraction diagram. 

8. Conclusion. 

As discussed in the preceeding section, the results of the refraction 
diagram agree well with that of the hand drawn method, so that the plots 
received from the computer can be used with as much confidence as those 
drawn manually. The decided advantage is that the computer product can 
be obtained many times faster. 

The coefficient of refraction and the wave height ratio as noted 
appear qualitatively correct and can be used with the reservation that 
the values have not yet been verified. 

9. Recommendations. 

Further development of programs similar to this will require a bet- 
ter representation of the bottom contours for more accurate results. The 
shallow-water depth values need to be very accurate for the proper deter- 
mination of all the parameters, since the depth determines the velocity 
of the wave at each point, and the other parameters are functions of the 
velocity. 

An accurate method within the approximations made is required for 


the calculation of the Beta parameter. The analytic approach of Kelvin's 


14 


method appears to be the more accurate method from the consideration of 
the necessary assumptions involved with the finite difference method. Іп 
the case involving steep gradients of the bottom contours, the finite dif- 
ference method as used here may prove inadequate. The time factor in com- 
pleting this paper prevented the successful completion of programming the 
Kelvin method as desired. 

Testing of the results is required to substantiate the procedure as 
a prediction method. Involved aerial photography is required to follow 
the wave crest refraction to the shoreline where the accurate prediction 
is most desired. Wave gages and other devices may be used to obtain the 
necessary data for the verification of the wave height values as predic- 
ted by the computer. A very intensive project will be required in order 


to verify the predicted values. 
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Figure 3. Computer drawn plot of wave crests 
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APPENDIX I 
Computer Program for Wave Refraction using Fortran 63 


on the C.D.C. 1604 Computer 


Program Title: WAVREFR 


Input Variables: 
АУ YV БЕ X and Y values of the first wave ray to be computed (feet). 


Bl, B2...... Values of the Beta coefficient (always equal to one in deep 
water). 


T .......... Wave period (seconds). 


DU κα. Angle measured from the direction of increasing X along the 
X-axis in the direction of travel of the wave crest (degrees). 


DOR аа 3 Number of rays that will be followed. 

ΡΕ Ες... Perpendicular distance between rays in deep water (feet). 

ΠΠ. ....:..: Time interval used for advancing the wave front (seconds). 
CRED ..... .. Distance between grid points of depth values (feet). 

ММ, ЫП ..... Number of grid points in the X and Y directions, respectively. 


DEP (I,J)... Depth values in the grid (fathoms). 


Output Variables: 


М... dum X and Y coordinates for a given wave crest and ray number 
(yards). 

GOREER ....« Coefficient of refraction for the wave at the point X,Y. 

ΠΕ Tele Ratio of wave height to the deep-water wave height at the 
point X,Y. 


NGO ........ Indicates that the ray has terminated (= 1), or is continuing 


(= 2). 
Variables in Common: (not previously defined) 
FK, FFK .... Values of ray curvature. 
CXY, CXX ... Values of wave velocity. 


PAR ee Constant = gT/2T. 
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BONES TE s uu Constant = 2T /T. 


MAX CC... V" Number of wave crest (= 1 for first crest). 
DEPTH TEL... Calculated depth at point X, Y. 

NORA ......; Number of ray being followed. 

O Е Indicates ratio of d/Lg. 


-1 for d/Lg > 0.5; 
=2 for d/Ig € 0.5. 
անա. ՐՈՆ Designates whether the last two curvature estimates for a 
given X, Y are less than .0001, whether the 13th and 15th 
estimates are less than .0001, or whether neither of the 
above is true (MIT = 1, 2, 3, for the respective cases). 
МІР ........ Used to determine the number of wave crest to be graphed. 
XXX, FX .... Values of X used in the graphing and printout. 
YYY, FY .... Similar to uses for XXX and FX. 
NNGO v, © The number of rays that have stopped. 
ХХ, ҮҮ, АА, А 
e oT Intermediate values of X, Y, and A. 
Summary of Program: 
The program reads the input variables, then the depth grid. PAR and 
BAR are computed for the wave period; A is converted to radians; and MAX 
and NORA are set equal to one. An initial value of CXY is found by test- 
ing the wave period. Control is transferred to the RAYN subroutine. When 
RAYN has determined the parameters at the first and second points along 
the ray, NORA is increased one, the XX and YY values computed for the next 
ray, and then RAYN is called again to compute the parameters. This pro- 
cedure is followed until NORA equals NOR, the number of rays. The graph 


is drawn for the first wave crest and the values printed. MAX is in- 


creased by one, and the procedure started over until all the rays have 
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stopped by either going off the grid, hitting the shoreline, or not having 
the values of the curvature converge. When all the rays have stopped NNGO 
is equal to NOR. The final graph is drawn which is the contour of the 
area utilizing the zero values of the depth grid and the program is сош- 


plete. 
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Subroutine Title: RAYN 
Variables of Subroutine: 


COREFA ..... Intermediate value of COREFR. 


ΕΠ  — ee Intermediate value of HHO. 
BABAR ...... Curvature used to obtain DEL A. 
NOGO жы е м, Storage value of NGO. 

ДАА τα... Storage value of A. 


Summary of Subroutine: 

RAYN calls DEPTFUN to obtain CXY, PDPX, PDPY, PDDPXY, PDDPXX, and 
PDDPYY. NGO is set equal to two if DEPTH is not zero, otherwise NGO is 
set equal to one. RAYN tests NGO of the ray for the MAX calculation to 
determine whether the ray will continue or not. If NGO = 2, KFUNCT is 
called to obtain the value of the curvature, FK. If NGO = 1, control is 
returned to WAVREFR. If the ray is continuing RAYN calls MOVE to project 
the ray to its next position. If the ray is not stopped in the MOVE sub- 
routine, RAYN next calls BETA to calculate the values of COREFR and HHO. 
As a final step RAYN stores the values of FK, COREFR, HHO, NGO, X, Y, 


PDPC, and AAA for use when the ray is again projected. 
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Subroutine Title:  DEPTFUN 
Variables of Subroutine: 


KER ........ Indicates errors in the solution of the simultaneous equa- 
tions for the quadric surface: 


- 1 no errors indicated in the solution; 


= 2 indicates that the matrix of values being solved is singular 
or nearly singular. 


МОК... Deep-water wave length. 
ШОУ an Ratio of DEPTH to ALO. 
BDB, BDEY EDDEXK, PDDPXY, PDDPYY 


..... First and second partial derivatives of the depth at X, Y 
with respect to the X and Y directions. 


Summary of Subroutine: 

The subroutine first determines the values of the closest grid point 
by truncating the values of XX/GRID + 1.5 and YY/GRID + 1.5 to give the 
correct value of the grid point. The quadric surface is fitted to the 
nine values of depth surrounding this calculated point, using the least 
Squares method, DEPTH is found by evaluating the quadric equation at X 
and Y values of the point. If DEPTH is positive, the wave velocity, CXY, 
is solved for by an iteration procedure of the equation described in the 
test using the principles of Appendix IV. The various partial derivatives 
are computed by evaluating the equation described in Appendix III. The 
method for the solution and the program for solving the six simultaneous 
equations resulting from the least squares method was written by C. B. 


Bailey and Mary Haynes of the USNPGS computer center programming staff. 
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Subroutine Title: KFUNCT 
Variables of Subroutine: 
PDPC, PDDPCC 


сый жаШ se First and second partial derivatives of the depth with 
respect to the wave velocity. 


PCPX, PCPY .. First derivatives of the wave velocity with respect to the 
X and Y directions calculated from relations in Appendix 
ІТІ. 
ЕК .......... Gurvature of the ray at the point X,Y. 
Summary of Subroutine: 
The subroutine calculates the values of PCPX and PCPY, and the curva- 
ture at the point X, Y from the equations of Appendix III and the text, 


respectively. The curvature is calculated only if LO = 2, where the wave 


1s in shallow water. 
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Subroutine Title: MOVE 


Variables of Subroutine: 


ΠΕΝ... Storage value of FKBAR. 

DD... Increment of distance advanced (DEL D - T CXY). 
ID Ee heres m XX - X. 

DIE NET аа YY - Y. 

DELA 0 АА - А. 

ABAR 39... (A + AA)/2. 


Summary of Subroutine: 

MOVE determines the X and Y values of the next point on the wave ray 
by an iteration process involving the curvature, the incremented distance, 
and the angle of the ray. The iteration is continued until the curvature 
estimates vary no more than .0001 from one another. Then the values of 
XX, YY, AA, and FKK are accepted for the new point. If the difference is 
greater, FKKP is set equal to FKBAR and the current FKBAR is used to obtain 
another set of values. If the iteration process stops before 15 itera- 
tions, MIT = 1. If the cycle stops at 15 iterations, and the difference 
between FKBAR and FKKP is less than .0001, MIT - 2, and FKBAR is defined 
as (FKBAR + FKKP)/2 for the last determination of XX, YY, AA, and FKK 
values. If the difference is greater than .0001 after 15 iterations, 

MIT = 3, and the ray is stopped. Control is sent back to RAYN. Inside 
the iteration loop DEPTFUN determines the parameters of CXY, PDPX, and 
the other derivatives, and tests the values of X and Y to determine if 


they are still on the grid. 
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Subroutine Title: BETA 
Variables of Subroutine: 
PGGPXX, PCCPXY, PCCPYY 


DP Vv Second partial derivatives of the wave velocity with respect 
to the X and Y directions. 


HSHOL ...... Shoaling factor. 
CGO Жее ыр Ratio of the wave velocity to the deep-water wave velocity. 


DD ......»... lncrement of distance along the wave ray between points n 
and n+l. 


B(i,j) ..... Values of the Beta parameter at the three points. 


Summary of Subroutine: 

BETA calculates the values of the coefficient of refraction and the 
wave height ratio at each point along the ray. The equations are des- 
cribed in the text. If LO = 1, BETA sets COREFR and HHO equal to one, 
and the Beta value equal to the previous Beta value on the ray. The 
shoaling factor is calculated from an equation determined from a poly- 


nomial fit to a curve of с/с, < 1 as described in the text. 
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PROGRAM WAVREFR 


WAVFONON 
WAVE REFRACTION PROGRAM FOR THE COMPUTATION OF CREST AND RAY WAVENO1O 
REFRACTION IN SHAOLING WATER BY LTe De STOUPPF, OCTOBER 1966. WAVE^020 
THE PROGRAM IS WRITTEN IN FORTRAN 63 FOR THE CDC 1604 COMPUTER МАУЕ0020 
SYSTEM. INPUT OF THE BINARY IS REQUIRED FOR THE DRAW SUBROUTINE WAVENNAN 
ON THE PRESENT SYSTEM AT THE US NAVAL POSTGRADUATE SCHOOL. WAVENNSO 
NOTE TO USERS ALL INTEGER INPUT ON THE DATA CARDS MUST BF WAVENNEN 

La RIGHT ADJUSTED FOR CORRECT RESULTS. | WAVENN79 
USERS MAY ENTER THE DESIRED NAMES FOR THE GRAPH TITLES IN THE МАУЕЛОВЯ 
ITITLE STATEMENTS AS DESCRIBED IN THE WRITEUP ON ТНҒ DRAW WAVENNSN : 

° SUBROUTINE WAVEN10N 
DIMENSION ITITLE(12) »XXX(100) sYYY(100) 9FX(100) sFY(100) Ὁ, 7МАУЕ0110 
COMMON/BLK1/X(100)»Y(100)»CORFFR(100)»HH0(100)>B(3» LOCI aia CUA МАУЕ0120 
1NOGO( 100) »sNGO»sAAA1 100) МАУЕ0120 
COMMON /BLK2/T y AsCXY sPAR, BAR, TIMES GRIDsSFKSMAX »;NORSMM, NN, DEPTH HHA; WAVE^140 
LCOREFAsNORA sLOsDIST WAVE9150 
COMMON/BLK3/DEP (1009100) »PDPX,PDPY, PCPX »PCPY »PDDPXX » PDDPYY WAVE^160 | 

| 13 PDDPXY ,PDDPCC,PDPC "WAVEO170 ` 
' READ2»XV»YV5B1»B2» T» ALS NORSDISTSTIMF,GRIDS,MM,NN WAVEN180 
2 FORMAT (2Ғ10.4192Ғ3.0>Ғ84.0»Ғ5.0,184»Ғ10.2»Ғ5.1»Ғ10.1»218) МАУЕ0190 
РВІМТ 19 WAVED20N 
19 FORMAT (1H195X92HXV 9 7X 9 2HYV 94X 5 PHRI 9 LX 92HR2 9 AX og LHT 9 2X a PHAL 9 2X 9 WAVE0O210 
13HNOR»3Xs4HDISTS4X5 4HTIMES 3X s 4HGRIDS, 3X, 2HMM$, 2X »2HNN/) WAVEN 220 — 
PRINT 2»XV5»YV5,B15»B25» T3 A1, NOR,DISTS, TIME, GRIDS,MM, NN ^ WAVFEn220 
READ 19t1(DEP(ISJ)»s Js 1»3NN) » I21, MM) WAVED 249 
1 FORMAT (14F5,0) WAVEN250 
^ 0022 1жізММ % 0022 J=lsNN МАУЕ 0260 

22 DEP(IsSJ) DEP(IS J) *6. WAVE^270 

, . .МАХя1 % А1х41%,017452329 $ PAR=32.2*T/6.283185 $ BAR=6.283185/T WAVEO2RO 

ρα 7 WAVEn29n 

. — LABELs4H $ DO 18 М=1,12 WAVEN3NN 

18 ITITLE(M)=8H τ 5 ITITLE(1)=8H STOUPPE $ ITITLE(4)=BHWAVE REFWAVED231N 

{ ITITLE(S)=BHRACTION $ ITITLE(6)=BHPROGRAM waAyrna2n 
ITITULE(8) 8HANGLE = $ ITITLE(10)=8H PFRIND WAVED 330 
ITITLE(9)=8H-35 DEG, $ ITITLE(11)=8H= 20 SFC WAVED 349 

~ -8(1»1)xB1 $ B(2,1):B2 $ XX-XV $ YY=YV $ NORA=1 $ A=A1 WAVEn350 
IF(T-10.)15»16,17 WAVEN 360 

. 15 CXY* 304,0 $ GO TO 3 WAVEO37^ 

-16 СХҮз 50.0 $ GO TO 3 WAVEN380 

. 17 CXY*80,0 $ GO TO 3 МАУЕ0290. 
2 6 В(1»ЫМОБА)-В1 % В(2»МОБА)-82 3 ΑΞΑΙ WAVEN4ON 
aces MVOEXVEDISTHSINFE (A) $ YV= YV+DISTCOSFIA)EXX=XV $ YY=YV $ 60 TO 3 VAVEN41N 
23 NORA=1 | WAVE^A420 
24 XXzX(NORA) $ YYsY(NORA) $ A=AAA(NORA) WAVEN430 
“° 8 CALt RAYN(XX»YY9NOR) WAVE^A44n 
CONTINUE $ NORA=NORA+1 $ SEE -NOR) 19,1094 WAVER ASN 

10 IF(MAX-110555,2^ | WAYEN4EN 
сезі! ҰЕЖМАХ-1)20%920»25 | WAVEN GTN 
25 IF(NIP-4)215,265,26 | WAVEOC4RN 
20 CALL DRAW(NOR»X»Y»1»2» LABELS ITITLE»s10000,5,100004,,0,0,2525999 sls WAVENA9N 
EAST $-GO TO 27 МАМЕ БОО 
26 CALL DRAW(NOR»X» Y s 2y2yLABEL I ITITLE 10000, ТОПОТ УЛОО 1 WAVEASIA 
1LAST) WAVEA^52^ 
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ei 


8 
7 


© 


12 
1H 
T3 
14 


30 
DT 


SUBROUTINE QAYNIXX. YY NOP ) 


COM 


O NZBLKLZFEFECCTUICOISDLLS 


NIE κι 
INOGO(190) »NGO> 


AAA(1CS 


Шао тти 


SA 1110911598111 
lS LNGO-?2 $ GC T2 158 
IMAN GO = NOGO NGO?) $ 52 T9 
LOOT GALL DEFPTFUN(XX YY) 

Ger 10 (194, πο ος 
104 NGO= D 090027203 
ШЕЙІРІМАХ- 13105,125,1^7 
WO Լ <FUNCT(AeFK) Հ ՞Չ 
ԱՏ K=FFY (NORE ) 
106 er. L E VY Y то 
ШЕ | ET ACK Xs YY) 
ЖМК ТИІ τ 027955 


т 


RETURN 
mE TURN 
-- perm 
E A 
РЕ AR AS Т 
MAS Тұ ТӨГІП 


ші rs 


ρωσ 
N 


]=HH4 t: wis 


аи wee 


- 
- 


ΙΙ τον 


Оа ae) MS 


Ve 


Mies c 


x ( 


Оор IT յին on EY 


ОРА) 


DI = X V 
. .. a 


a CM HS NS RIED БАРЕТ,ІТІТЕНМ100272..12020.Э 020, A EN SA 
IE ST) Բ MAVFA54^ 
mni MAVEARRA 
ПИЕ МІР+1 5 2 τ. ινεύλχ-] WAVE NE KN 
268932 у-1.409 ՀՐԱ ЕЕ сива, MAVEA^RTA 
CONTINUE WAYENS2N 
արել 28214 . WAYEN.SON 
FORMAT(//17HNUMBER OF CREST = 14/) YIAVF ^ 6AN 
PRINT 6 WAVEN619 
FORMAT (/4X91HX>)12X»1HY98X9€=ACOREFRs33X93HHHNs 7x9 3HNGO+10X91HX312X M"AVE^62n^ 
11НҮ 9 8X s GHCOREF 25 3X 9 SHHHO 9 7X9 BUNGOS) WAVEN630 | 
DON J-J s NORIS | MAVERAAG^ | 
ΓΝ B.EXCULTSEYTCJYSCOSEES (LI) ο ος η Εκ J:l1 y, Ey ( SIE E δρ B 
ICOREF^?(J41)5HH2(J41) , "050 (J*1 ) WAVEAA6^ 
ИПГАТ(2(Е10.227Э3УХ>ҒІС,2,2Х,Ғ55,?,?У,СЕС,?,3У,1297Х)/) WAYENKTA 
CONTINUE WAYENEAN | 
2 І-1,М0? % DO 9 J=252 WAYENGOR | 
Επι ισα ντ) = ММО Е Բյ: κε  ΝΟΡ WAYENTAN | 
IF(NOGO(K)-1)12512511 B WAVEA^7]^ 
NNGO=NNGO+1 WAVE^T7T2n | 
CONTINUE $ IFLHNSO-NOR) 13,14, 14 է WAYENTAN ' 
το ο” ' WAVEn74n Ë 
CONTINUS WAYERNTEN | 
ΜμΞ1 8 ρο31]-15»»'! 8 DO31Js1,NN « IF(DEP(ISJ))315»30,31 Hau σε E 
XXX(M)=(I-1)*¥GRID © YYY(“)=(J-1)¥*GPID $ M=M+1 $ N=M-] WAYENTTN | 
CONTINUE WAYEA TAN 
ӘРІПІ ГЗАЛШИБХХХ»ҮҮУ»30.51Ла8ас|,ТТТІТІ 5,1000, 1000.47, 2 es ος 
ІШ БАСТ) WAVER AAA 
END WAYENL]IN '| 
| 


ΏΛΥΝΛΑ2Λ 
RAVNARAN 
RPAYNQ 849 


RrE5(100) ,HP^(1^22) 9P(2, 10^) ,FEK(10^), 


ABC MMI NNI DEPTH HHA, 


COMMON/BLK2/T»A»sCXYsPAR,BAR,TTME,CPID,FK,.MAYX > Աո eG of 

NSCHSECA.NORAsSLOsDIST ΟΑΥΝΩΒό6Ω | 
COMMON/BLK3/DEP(102,100) S PDPX,PDPY ,PCPX,PCDY,PDDPXX, PDDPYY ΟΛΥΝΛΑΤΩ | 
1: PDDPxvY ,PDDPCC,^DPC | олумладл | 


RAYNNA9nN . 
RAYNNQNA | 
DAYNNQ1N | 
DAvM^O2n Pl 
O AMAS 20 
DAYN^OA^ 
PAVNAQEA j 
DAYNNOAr 
RAYN^OT7^ |. 
ΘΔΥΝΩΟΑΘΩ 
QAVNNADAN | 
олум] тил | 
DAvM1^1^] 
O^vYN1^2^. 


.o. C EEWEOTA 4 HHWN-OVO t col To 163 


1 
ж 


<> 


6 


4414 A ( KUNDEN ) 


“ ү(мбоку-суу FT PPATRNIADAYHA 


NAYRI AA | 
navuının | 


ы гут 1 с ^Ш 





30 


աե ԼԶԸ 2(6» os ο) DUC DE I) 

BE cixi/xclo OysY(12C) COR CER ( 1007) μα ο δις ο το ο ο. 
1М960(1003»%60»ААА( 1220 | | 
COMMON /BLK2/T sA a OXY s PAR BAP TIYE GRIDI FLI MAX NOR Mi NN SDED TH gHHA g 
B E ANORA LO DIST 
COMMCN/BLK3/0EP (19000109) . РОБУ гпрү,эсрх „РсРрү s,PDDPXY ,PODPVY 
1»PDDPXY»PDDPCC,PDPC 

Meee eng DLS 4 /FF KK (100) Ü CL 2 6 C(—(6) "IT DSL V DEC Y 

K ΙΙ ç K=XR+1,5 {f YR=aYY/TOID € [=VD+),F 

տն 2-120)2(0(/4-1)-Հ) )5227»7 


17-12.) ((ԳՈ-1)-Լ))52:65»6 
6 DOSJ=A1,6 Ֆ 0091-5196 
9 D(I>»J)=0.06 $ 008 І=1,6 
и ЕТТ)-02.0 
M=L=1 $ MA=L+1 S M=K-1 5 M n 
0010 J=N,MA Ֆ 0Օ1Շ 15 =P a DLA “O( 1)=€1-1)4GR1ID Ρις уе = DD 


D(1s1)=9e ο ο ociosa οι”. 17+2( ЈО (АУТ) ЕСО ТО 
ο pUo5s1)-D(5.1)4TOC(I)**2$05(6,1)23D(5,1) *P( J| **2$ FE ( T1) ECT) FDEP 
ENERO LUIS 2)sD(2»1I$0(252)2D2(5,120€D(352)2D( 4,9 1)$D(452)50904.92) OCT 
Ա Թ J) S$D(532)=0(552)+0(1)**#29D(552)=0(652)4+90(1) *XP(J)*X2$F(2)s 
DEIN UP)A*DEP(LIS$999$5(153)20(3,120$0(25,3)2^(4541) 80( 3, 23) SD (65,42. SP( 43 
ο DUAS 3) OCT) *XPUJ)**2$D( 5,3)-^5(4,42)$D(6.$3)2D(65 3) P CJ) 5X3 $F Ἅ)}Ξ 
DEEP *DEPUISJ)SDC,I1.S4)2D(4,1)$0(2,4)20D(45,72) 5DUC 2.44 τος τοις, 


74)=D(494)+0(1)**2*P (J)**25D[(594)=D(594)+0(1)**3x*D(J)$D(694)=D(64 


8)4O0CI) *PCJ) **3$E(4)2E(4) OCT ) DC J) X5DEP(1,0)$^(1,5)20(5,1)$D(25 5) 
9D( 5»21$D( 3»5)2D(4»2) $5D( 4»55)20D(65,4)$D(5,55)20(5,5 5) OCT ) **4$0( 5,5) 
ο oe (5) =E(5)4+0(1)**2*DEP(] -J)$5D0(196)=D(631)5D(2245)=D6s2) 503 
256) =0(6 93) 0р (496) = 0 (65 4) $0 (5,6) =0 (4 4)$9(0S957=D0696) +P U2) АН 
3=E(6)+P(J)**2*DEP(I9J) | 

ШО CONTINUE Р 
τ. 50 344-21.»6 5 KP=0 Ֆ Z=0.0 $ DOLZN=M56 5 TC(7 TDN 
ШІ! 11512,12 

11 2-АВ5Ғ(О(М»М)) $ KPzN 

W2 CONTINUE $ IF(M-KP)13920920 

թի Ն ի չիթ ի %7/-0(М».)) $ D(Ms3J)=D(KP9J) 

14 О(КР» 2) =2 

КОТЕ САВБЕСОСМ, М) ) 00001 ) 50, 50, Зп 

30 ΙΕ(Μ-6) 2194040 

31 LP1=M+1 $ DOZB4N=LP155 5 IF(D(N SM) ) 32934532 

BZZRATTO=EDIEN»MY/DIM>M)S DO?3J=LP1,NPM 


B3 DNsJ)zDUNS J) -RATIO*D(Ms$,J) 

34 CONTINUE 

աւ. 151:6 $ IIs7-I 5 ре Ф 5-С.0 3 TEC(TISS 1415055 
143 


4] IIPl*II4*1 $ DO 42N=IT 19155 

42 SESEDOIISNO*C(N) 

ES C(II) SUDINLETS,JPN)-S)/D(CIISII) TF KESSI $30 p 7 

50 KER=2 | 
PRINT 54» ХХ»ЭУУЭКЕК 

ЕШ FORMAT(/3X95HXX = 1PE20.8»3Xs5HYY = 1PF27.2, 4X» S5HKEP = T2, 
AT MATRIX SINGULAR/) 

տն նց (52153)KER 

53 RETURN 


on 


DEDTl1^6^ 
AED TAZA 
DEPT1^80 
DEoT1^90 
DEPT11^^ 
DEPTI11N 
Ծ657112Դ 
ԹՔԾ71112ԴՊ 
DFDT114^ 
NFPT115N 
DEE IIS 
Dep τη 
DEPT118N 
DEPT119^ 
DEPT120^ 
БЕРТ τρ 
ШЕРДІ ОВ 
DEOTI 239 
DCcPT124^ 
ηεοτ12τω 
Բրա 22A 
DEDT1279 
ΏςΩΤΙ 22 
DEPT1290 
DEPT! 200 
ΠΕΡΤΙ341Ω 
DEPT1320 
ncpT123n 
DEPT1340 
DETI 250 
DFPT12350. 
DFPT1279 
NFOT1 294 
DcPT139^ 
перт] 407 
NFPT1419 
DEDT1I42N 
DEPT1439 
DEPT1440 
DEPT1450 
DEPT146N 
DEPT1479 
ΠΕΡΤΙ4ΑΩ 
DEP ОП) 
DEB MSC 
DEPTISIN 
DcPT1520 
DEPT lean 
DEPT154^ 
DEPT1550 
DEPT156N 
DEPT1570 
DEPT158^ 


өз 
No 


62 
64 
81 


E 
6 


1 


CONTINUE 

DEPITE (])+C(2) СТ) УУСС) ОИ 
IF (DEPTH) 62562565 
ALO=PAR*T 5 DLO=DEPTH/ALO 
ӘЛЕГІ $ LO=] “GeNiGO=2? 3 FKk=0 D OMM 
DO60M=1»5C $ CXX=PAR*TANHF (BAP»*DEPTE/CXY) 
IF(ABSF(CXX-CXY)=0.01)61>61960 | 
CXY2(CXX-CXY)?t, 5 

CONTINUE 3 NGO=2 $ LO=2 
PDPX=C(2)+C(4)*YY+2*C(5)%*XX 
PDPY=C(3)+C(4)*XX42*C(6)*YY 
PDDPXY=C(4) $ PDDPXX-2*C(5) 
GO TO 64 | 
NGO=1 $ 60 TO 81 
CONTINUE 

RETURN 

END 

SUBROUTINE KFUNCT (AsFK) 
COMMON/BLE1/Xt100) »Y C160Y s CORFFR(10015HHC (1003 38(3, 100) SEFK (10^) » 
1М060(100) sNGOsAAA( 100) 

COMMON/BLK2/T UUs CXY sPAR s BAR sTIME sGRIDs VV oMAX sNOR 9 M4» NN a DEPTH HHA 
1»COREFA»NORAsLOsDIST 


Merete PEG) EY YER 


% ІБ(СІ2-0.5)52,623,66 


$ PODPYY=2*C(6) 


τν O S DUK3/DEP”(100 I 100)sPDPX POPY s PCP Q C PY sPODDXXyPDDPYY 

Jy Pi C PDDPCCyPDPC 

Ὅο ο ο δτι ο | E 
R12CXY/32.2 $ R2-R]**2*BAR**2.$ R3-R]**5*BARGX4 Տ Y =RIAATARARAAS 


PDPCz2,.*R1-4,.*892/3,.--64.*R3/5,^*AR,*R4/7, e 
PDDPCCzZ(2,.*R14T4,.*X82«46, *82*8,*X94) /CXY 
PCPX-PDPX/PDPC $ PCPY-PDPY/PDnrc 

ԱՆ Թ  SINPT (A)=PGPY*COSF(A))/CXY $ 
FK=08 

RETURN 

END 

SUBROUTINE MOVE (X Y) 
COMMON/BLK1/U(100) »V(100) » COREFR 
1NOGO(100)»NGO»>»AAA(100) 

COMMON/ BLK2/T sAsCXY sPARsBARs TIME sG 
1COREFA,NORA,LC»DIST 
COMMON/BLK3/DEP(10C$,100) » PDPX,PDP Y, PCPX,PCPY ,DPDDDPXX , DhDPYY 
1»PDDPXY»PDDPCC,PDPC 

COMO pPEK4&/FFKK1 (100) + DLU0 C (6) s IT | DSL X,DEL Y 
FKBAR=FFKK(NCORA) 

IF (MAX-1) 19194 

FKBAR=FK 


СОТ -4 


(1201: H501106)9B.3» 100 ΤΕ κ», 


RIDsF Ko MAX з. МОВ о ММ ao NN a DEPTH gHHA, 


4 MIT=1 


22 


21 
po 


DEL DsTIME*CXY 
Come (22521)L0 
XX=X+DELD*COSF(A) 
GO TO 6 

DONC CTTsS 15 

DEL A=FKBAR*xDEL D $ 
1 COSFCABAR) 
ՏՕ ТО 


P oYYSYTDEED*SINECA)OEEENEN 5 ԲՇԸՀՀՍՀ» $ FRAIS, 


AA=A+DEL A $ ABAR=A+e5*DEL A $ 
$ DEL Y=SFL D*¥SINE(ARAR) 
(101>61МІТ 


hae A ο 
Տա Σχ. κ S$ YYzY+NEL Y 


32 


EPIL Ρ 
DebPTTen^ 
ОЕР ТИХО 
ОЕР Ա 


DEPT162^ 


DEPT164^ 
DFPT165N 
DEPT166N 
DEPT1670 
DEPT1680 
DEFT1690 
DFPT1700 
D: 3271 10 
DEPT172^ 
DEPT173^ 
DEPT TAN 
DEPT1 75^ 
KFUN1 760 


KFUN177^ 


KFUNT 789 
KFUN1 799 


KFUN1800 © 


KFUN1 819 


-KFUN1822 


KFUN1830 
K FUN1840 
KFUN185N 
«ԲՍՎ1860 
KFUN187n 
KFUN1880 
KFUN1 890 
KFUN19۵0 
KFUN1910 
MOVE1920 
MOVE1930 
М0УЕ 1940 
MOVE1959 
MOVE 1960 
MOVE} 979 
MOVE198N 
М0У51990 
MOVE29N4 
MOVE 2910 
MOVE 2920 
MOVE2N3N 
MOVE 2049 
MOVE 2050 
MOVE 260 
MOVE 2070 
MOVE 2089 
MOVE 2N9N 
MAVE2 TAN 
"OV£211^ 


- 





БАР МЕРТЕШ ХХ»ҮҮ) 
CALL KFUNCT(AA,FKK ) 


IFCIT-1325,37,9 
FKKPP=FKBAR 

IF (MAX=-1) 79799 
IF(IT-1)20+2099 
IF(ABSFIFKKP-FKBAR)-4.,00001)656,20 


FKKP=FKBAR 
IFLABSFIFKKPP- ҒКВАВ)-.00001718, 18,17 


$ GO ТОСЗБИТОТМЕВ 
$ FKBAR=e5*(FK+FKK) 


MIT=3 $ NGO=1 $ GO TO38 


FKBAR=.5* (FKBAR+FKKP) 
NGO=2 $ GO 


ro. 8 


Х=ХХ $ Y=YY $ AZAA $ FK=FKK 


CONTINUE 


FFKK (NORA)=FKBAR 


RETURN 
END 


-* SUBROUTINE ВЕТА(ХХ»ҮҮ) 


7 


COMMON/BLK1/X(100)»Y(100)»CORFFR(100)+HH0(100)»B(3» 100)»FFK(100)» 
Τι «ՅԱ AAA( 100) 


S.MITs24$ GO ТО 19 


MOME 2120 
MOVE2130 


MOVF2140 


МОУЕЙМФО 
ΜΌΝΕ ΕΟ 
MOVE2170 
MOVE2180 
MOVE2.190 


MOVE220N 


МОУЕ2210 
М0УЕ2220 


MOVEZ220 


MOVE224^ 
MOVE2250 
MOVE2260 
MOVE2270 
MOVE2 280 


BETA2290 ` 


BETA23N0 
BETA2310 


° COMMON/BLK2/T sAsCXYsPARs ВАК» ТІМЕЭСКІПОЭЕК»ЭМАХ»МОБ»ММ»ММ»ПЕРТНЭННА» ВЕТА2220 


“ACOREFAsNORAsLO> DIST 


COMMON/BLK3/DEP(100+100)»PDPX+,PDPY » PCPX, PCPY »PDDPXX sPDDPYY 


1 sPDOPXY sPDDPCC »PDPC 


COMMON/BLK4/FFKK(100)5DLOsC(6), MIT. DEL Xs,DEL Y 


GO TO 


GO TO 7 


(5»6)Լ0 
COREFA=le 


$ HHA=le 


$ B(3>NORA)=8(2»NORA) 


PCCPXX=PDDPXX/(PDPC+PDDPCC) 
FOCPTYSPDDPYY/(tPDPC-*PDDPCC) 


. RCCPXY=PDDPXY/LPDPC+PDDPCC) 
Pz(-COSF(AY)*PCPX-SINF (A) *PCPY)/CXY 
)SINF(A)*COSF(A)*PCCPXY«*COSF(Ay*2*PCCPYY)/CXY 


DD=SORTF( (DEL X)**2+(DEL Y)**2) 
ВЗ МОКА ) = В(1 + МОВА ) + (Ржрр- 2, ) +В(2 5 ЧОЗА ) # (40-2, #00260) ) 


1/t2.+P*DD) 


СОКЕҒАғ1./<ОКТҒ(АВ5Ғ(В(2»МОКА))) 


ΗΠ 3 HSHOL=3.2519-12.8150*CC0+28.8112*CC0*2= 29.92 5 70ECO 
1**3+11,6815*CC0**4 $ HHA= COREFA * HSHOL 


RETURN 
END 


335 


$ 


Ф СОРЕЕВг1. 


Փ HHB=1. 


Qs(SINF(AY)**2*PCCPXX-2.* 


BETA2330 
BETA2340 
ВЕТА2250 


BETA2360 ` 


BETA2370 
BETA2380 
BETA2390 
BETA2400 
BETA2410 
BETA2420 
BETA2430 
BETA2440 
BETA2450 
BETA2460 


BETA2470 © 


BETA2480 


ВЕТА2490. 


ВЕРА2600 
ВЕТА2510 
ВЕТА2520 


ται 


. 
"e 


253 


Example of Input 


X Y 
4500. 33000. 


82 Т А 
1. 20.-35. 


81 
le 


NOR 
18 


DIST 
1500. 


TIME 
30. 


GRID 
1500. 


MM 
52 


NN 
56 


440, 
340. 
550. 
510. 
360. 
290. 
540. 
610. 
295. 
240, 
560, 
540. 


460. 
345. 
490. 
460. 
405, 
297. 
500. 
560. 
370, 
250. 
540. 
590. 


215» 
350. 
470. 
400. 
450. 
300. 
515. 
515. 
370. 
250. 
550. 
570. 


505, 
365, 
460. 
315. 
435. 
305. 
515. 
200. 
325, 
265, 
560. 
430. 


510. 
375. 
415, 
275. 
445, 
315% 
4800 
315. 
325. 
275. 
550. 
300. 


470. 
390. 
380. 
300. 
380. 
360. 
4456 
275. 
307. 
305. 
540. 
310. 


425. 
455, 
3356 
220. 
330. 
420. 
440, 
234. 
260. 
405-6 
550. 
280. 


355, 
560. 
400, 
140. 
270. 
490. 
500, 
190, 
190. 
470. 
605. 
260. 


300. 
700. 
460. 

80. 
235. 
6556 
540. 
110. 
170. 
570. 
620. 
170. 


330. 
850. 
520. 

55. 
250. 
810. 
592. 

63. 
175. 
750. 
640. 

Ը. 


370. 
950. 
560. 

53. 
280. 
930. 
630. 

53. 
190. 
900. 
610. 

55. 


360. 
865. 
595. 

51. 
280. 
870. 
630. 

50. 
190. 
885. 
570. 


340. 
705. 
620. 

51. 
260. 
750. 
600. 

50. 
180. 
800. 
530. 


330. 
630. 
590 

51. 
2706 
620. 
590. 

50. 
205. 
630. 
510. 


240. 
190. 
610. 
4800 
150. 
162. 
660. 
450. 
80, 
110. 
760. 
440. 
61. 
80. 
770. 
4256 
blo 
630 
560. 
420. 
38. 
60. 
450. 
430. 
31. 
59, 
360. 
500. 
30. 
58. 
2400 
460. 
13. 
56. 
155. 
220. 
18. 
55, 


250. 
195, 
590. 
520. 
195. 
170. 
650. 
490, 
80. 
120. 
720. 
410. 
60. 
85. 
750. 
490, 
47. 
656 
610, 
510, 
38. 
62. 
465, 
520» 
320 
60e 
380. 
470. 
28. 
59, 
280. 
410, 
24. 
57. 
180. 
370. 
20. 
566 


2356 
210. 
590. 
560. 
200. 
180. 
6406 
550. 
90. 
140. 
715. 
550. 
60. 
100. 
735. 
550. 
48. 
70. 
590. 
515. 
38. 
63. 
490. 
870. 
33. 
61. 


195, 
230. 
615, 
515. 
110. 
205. 
660. 
525, 
75. 
160. 
725. 
490, 
58. 
120. 
715. 
500. 
49, 
7. 
600. 
500, 
41. 
67. 
510. 
430, 
38. 
62. 


420.415. 
410,360, 


33. 
60. 
325. 
300. 
26. 
58. 
2400 
300. 
240 
58. 


36. 
60. 
225. 
260. 
28. 
59, 
2400 
2400 
250 
580 


2256 
250. 
600. 
405. 
125. 
210. 
660. 
450. 
73. 
1756 
730. 
415. 
60. 
130. 
665. 
4456 
52. 
82. 
570. 
465. 
46. 
70. 
470. 
415. 
40. 
63. 
400. 
350. 
35. 
62. 
325, 
260. 
28. 
59, 
230. 
240. 
25. 
57. 


200. 
280. 
610. 
370. 
120. 
210. 
670. 
380. 
70. 
195. 
710. 
370. 
60. 
120. 
635. 
375. 
54. 
85, 
5400 
410. 
48. 
70% 
465. 
820. 
43. 
64, 
390. 
370. 
356 
62. 
330. 
295. 
31, 
60, 
2306 
2100 
260 
59. 


190. 
350. 
620. 
320. 
115. 
290. 
690. 
325. 
72. 
240, 
680. 
325, 
61. 
160. 
620. 
315. 
53. 
110. 
525. 
283. 
51. 
80. 
460. 
240, 
48. 
65. 
390. 
380. 
38. 
63. 
340. 
300. 
34. 
60. 
260. 
180. 
27. 
60. 


160. 
427. 
660. 
260. 
95, 
380, 
685, 
270. 
70. 
290. 
6556 
285. 
63» 
230. 
610. 
295, 
She 
150. 
520. 
265, 
53, 
95, 
460. 
270. 
50, 
70. 
410. 
375. 
436 
646 
3306 
290. 
28. 
62% 
270. 
200. 
330 
60. 


34 


130. 
490, 
650. 
190. 
80. 
4556 
650. 
215. 
70. 
390. 
620. 
235. 
60. 
280. 
570. 
265% 
56. 
185, 
500, 
235, 
54% 
115. 
4400 
254, 
Sle 
80. 
380. 
340. 
46. 
65. 
310. 
300. 
42. 
64% 
270. 
2000 
360 
62 


140. 
670. 
620. 
80. 
98. 
580. 
585. 
120. 
70. 
475. 
5400 
1706 
60. 
380. 
485. 
190. 
59, 
250, 
460, 
210, 
56° 
160. 
415. 
260. 
54. 
86, 
365, 
310. 
49, 
73. 
290. 
300. 
45, 
66, 
280, 
190, 
40. 
63. 


160. 
820. 
560. 
60. 
105. 
740, 
510. 
656 
12. 
610. 
480. 
85. 
656 
460. 
440. 
110. 
60. 
340. 
410, 
145.6 
51. 
210, 
395, 
250. 
55, 
140, 
350, 
315, 
Sle 
80. 
305. 
290. 
4B, 
70. 
325, 
210. 
44. 
63. 


50. 
170. 
860. 
505. 

49. 
110. 
785. 
480. 

50. 
80. 
670. 
445 
65. 
65. 
540. 
415, 
80, 
60. 
395, 
380, 
105, 

57 
300, 
385, 
215, 

56. 
185. 
350. 
300. 

53. 
120. 
320. 
310. 

51. 
120. 
370. 
220. 

49, 

85. 


50. 
170, 
820, 
480, 

49, 
115, 
850. 
465. 

50. 
85, 
750, 
430, 
50, 
75» 
630. 
385, 
65, 
60. 
465. 
360. 
60. 

59, 
370, 
3700 
160. 

57. 
285, 
350, 
195, 

Dos 
185, 
405, 
300, 

54 
160. 
440 0 
3000 

51. 
120, 


50. 
180. 
680. 
470% 
49. 
135, 
720. 
440. 
50. 
95. 
815, 
420. 
50. 
75, 
775. 
375. 
50. 
62. 
550. 
365, 
52. 
60. 
450, 
360, 
55. 
58. 
385, 
3520 
120. 
56. 
260. 
450, 
220. 
54.6 
160. 


470. 


260. 
53. 
105. 





Example of Output 


NOR 


> 
о 
о 


XV YV B) 82 T Al DIST TIME GRIO MM NN 
4500.0 33000.0 1 1 20 -35 18 1500.00 30.0 1500.0 52 56 
NUMBER OF CREST = 1 
X Y COREFR HHO NGO X Y COREFR HHO 
2339.60 10412.11 1.000 1.000 2 \ — 2626.38 10821.68 1.000 1.000 
2913.17 11231.26 1.000 1.000 2 3199.96 11640.84 1.000 1.000 
3486.75 12050.91 1.000 1.000 2 3773.54 12459.99 1.000 1.000 
4060.33 12869.56 1.000 1.000 2 4347.11 13279.14 — 1.000 1.000 
4633.90 13688.72 1.000 1.000 2 4920.69 14098.29 1.000 1.000 
5207.48 14507.87 1.000 1.000 2 5494.27 1կ917.ԿԿ 1.000 1.000 
5781.05 15327.02 1.000 1.000 2 6067.84 15736.60 1.000 1.000 
635% „63 16146.17 1.000 1.000 2 6641.42 16555.75 1.000 1.000 
6928.21 16965.33 1.000 1.000 2 7215.00 17374.90 1.000 1.000 
NUMBER OF CREST = 2 
x Y COREFR HHO NGO x Y COREFR HHO 
3179.19 9824.22 1.000 1.000 2 3465.98 10233.79 1.000 1.000 
-. 3752.77 10643.37 1.000 1.000 2 4039.56 11052.9% 1.000 1.000 
4326.35 11462.52 1.000 1.000 2 4613.13 11872.10 1.000 1.000 
4899.92 12281.67 1.000 1.000 2 5186.71 12691.25 1.000 1.000 
5473.50 13100.83 1.000 1.000 2 5760.29 13510.%0 1.000 1.000 
- 6087.07 13919.98 1.000 1.000 2 6333.86 14329.55 1.000 1.000 
6620.65 14739.13 1.000 1.000 2 6913.46 15160.19 1.000 .905 ` 
7182.28 15566065 1.000 .899 2 7454.38 15986.51 1.000 .897 
7728.66 16404.84 — 1.000 .895 2 8001.84 16823.95 1.000 .895 
NUMBER OF CREST = 3 
E Y COREFR HHO NGO X Y COREFR HHO 
8005.32 9245.76 1.000 .902 2 4292.80 9654.85 1.000 .899 
4577.71 10065.74 1.000 .897 2 4858.41 10479.58 1.000 .897 
5180.28 10890.72 1.000 .897 2 5423.42 11304.73 1.000 .898 
5711.14 11713.66 1.000 .898 2 6001.92 12096.33 1.000 .898 
6307.35 12513.11 1.000 .903 2 6587.17 12931.41 1.000 «903 
6854.59. 13354.55 1.000 .898 2. 7124.80 13775.74 1.000 .896 
7398.65 14197.17 1.000 .895 2 7694.29 14666.16 1.024 .917 
7930.90 15056.20 1.004% „899 2 8191.41 15480.48 1.032 .923 
8446.16 15904.45 1.03] .922 2 8727.76 16319.95 1.009 .903 
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FS һу һу Թ Թ Թ Թ ο N 


N N N N Թ Թ һу Թթ м 


т 
e 
o 


z 
e 
o 


N N N N N N км к 


NUMBER OF CREST = 22 


X 
10268. 28 
11438.10 


12920.26 . 


17781.43 
18531.25 
19101.90 
19789.68 
20588.90 
20812.04 


NUMBER OF CREST 


x 
10264.28 
11438.10 


77 12920.26 


18177.11 
18949.52 
19523.12 
20218.%% 
21014.87 


22 21212.27 


NUMBER OF CREST 


x 
10264.28 
11838.10 
12920.26 
18502.09 
19320.29 
19879.31 
20591.71 
21375,95 
21548.37 


Y 
4616.24 
4352.73 
2737.33 
2925.94 
4023.30 
4934.51 
6181.58 
7626.71 
8068.66 


23 


Y 
4616.24 
4352.73 
2737.33 
2611.34 
3760.33 
4669.33 
5951.79 
7416.37 
7844.86 


= 24 


Y 
4616.24 
4352.73 
2737.33 
2350.98 
3526.00 
4454.54 
5763.85 
7248.08 
7673.71 


COREFR 


0 

0 

0 
0467 
. 882 
«757 
1.139 
«930 
1.437 


COREFR 


0 

0 

0 
«38 
«882 
«725 
1.143 
. 942 
1.365 


COREFR 


0 

0 

0 
e415 
„863 
«708 
1.124 
. 945 
1.332 


HHO 


2486 
2924 
«791 
1.199 
„986 
1.542 


HHO 


.488 
.962 
.807 
1.270 
1.068 
1.584 


HHO 


.536 
1.084 
„915 
1.408 
1.397 
1.864 


с 
ο 
ο 


Ի) Ի) Ի 9 Ի» Իշ ա տ ա 


= 
с 
O. 
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2 
e 
o 


N NN NNN œ ~ ου 
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X 
1085 3. 02 
12072.87 
15666.05 
16546.11 
18968.07 
19511.47 
20449.56 
20819.12 
21047.31 


X 
10853.02 
12072.87 
15818.45 
16808.11 
19386.59 
19940. 13 
20883.05 
21227.05 
21454.05 


X 
10853.02 
12072.87 
15936.15 
17039.75 
197Կ3.52 
20297.47 
21240.91 
21569. 75 
21784.38 


Y 
4550.71 
3316. 63 
1290. 81 
1779. 73 
4701.04 
5662. 92 
7333.27 
8054.09 
8517.89 


Y 
4550.71 
3316.63 

860.81 
1362.17 
4438.11 
5417.85 
7120.80 
7848.73 
8311.97 


Y 
4550.71 
3316.63 

482. 73 

997.51 
4220.11 
5212.49 
6950.72 
7698.52 
8176.83 


COREFR 


0 

0 
«270 
8.942 
. 769 
.859 
1. 204 
2.473 
1.045 


COREFR 


0 
0 


2252: 


6.696 
«756 


„850 


1.131 
2.254 
1.031 


COREFR 


0 
0 

. 238 

5.649 
.753 

.843 

1.088 
2.148 
1.041 


HHO 
0 
0 
„291 
9.377 


2 805 


„900 
1.271 
2. 659 
1.125 


HHO 


«287 
7.344 
.840 
. 949 
1.285 
2.625 
1.228 


HHO 


„296 
7.103 

. 964 
1.131 
1.493 
3.043 
1.774 
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e 
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1 
2 
3 


SUBROUTINE DRAW (NUMPTS» X>» Y» 


A GENERAL CURVE 


PROGRAMMER 
DATE 
SYSTEM 
OUTPUT 
NOTE 


APPENDIX II 


INPUT ARGUMENTS --- 


le 


le 


3e 


4e 


De 


6» 


NUMPTS 


MODCURV 


ITYPE 


LABEL 


MODCURVs ITYPE, LABEL, 00000000 

ITITLEs EXSCALEs YSCALEs IXUP, IYRIGHT» 00000010 
МООЕХАХ» МООЕҮАХ» ІМІПОЕ» IHIGH, IGRID, 00000020 

LAST) 30 
0 40 

50 

DRAWING AND POINT PLOTTING SUBROUTINE 00000060 
Je Re WARD i 00000070 

РЕВ, 19645 REVISED JUNE 1965 00000080 

FORTRAN 60 00000090 

LOGICAL TAPE NUMBER 8 00000100 

ASTERISKS MARK CHANGES FOR FORTRAN 63 00000110 

000 120 

00000130 

000 140 


NUMBER OF POINTS TO BE PLOTTED. 


THIS MUST ALWAYSO0000150 
BE AT LEAST 2» AND MUST NOT EXCEED 30 FOR POINT 00000160 


PLOTTING» OR 900 FOR CURVE DRAWING. 00000170 
000 180 

ARRAY OF X-ORDINATES. DIMENSION AT LEAST EQUAL 00000190 
TO NUMPTS AND NOT MORE THAN 900 IN CALLING 00000 200 
PROGRAM, 00000210 
000 220 

ARRAY OF Y-ORDINATES. DIMENSION AS FOR THE 00000 230 
X ARRAY IN THE CALLING PROGRAM, 00000 240 
000 250 

CONTROLS THE NUMBER OF CURVES» AND/OR SETS OF 00000 260 
POINTS» ON EACH GRAPH, THE CODES ARE 00000270 
0 ONLY ONE PLOT ON THIS GRAPH 00000 280 


1 FIRST PLOT ON MULTI~PLOT GRAPH 00000290 

2 INTERMEDIATE PLOT ON MULTI-PLOT GRAPHO0000 300 

3 LAST PLOT ON MULTI-PLOT GRAPH, 00000310 

000 320 

CONTROLS THE TYPE OF PLOTe THE CODES ARE 00000330 
0 STRAIGHT LINES JOIN SUCCESSIVE POINTSO0000 340 
(STANDARD CURVE DRAWING) 00000 350 

1 POINTS PLOTTED WITH CROSS (X) 00000360 

2 POINTS PLOTTED WITH PLUS (+) 00000370 

3 POINTS PLOTTED WITH SQUARE 00000380 

4 POINTS PLOTTED WITH DIAMOND 00000390 

5 POINTS PLOTTED WITH TRIANGLE 00000400 

WHEN POINTS ARE BEING PLOTTED (ITYPE=1 THRU. 5)+00000410 
THE POINTS ARE NOT CONNECTED. 00000420 
000 430 


IF A CURVE IS BEING DRAWN (ITYPE = O)» LABEL IS 00000440 
A 4-CHARACTER BCD CURVE IDENTIFIER WHICH WILL BE00000450 
REPRODUCED AT THE END OF THE CURVE» FOR EXAMPLE,00000460 
LABEL = 4H ΟΝΕ ο IF POINTS ARE BEING PLOTTED, 00000470 
LABEL IS АМ 8-CHARACTER IDENTIFIER. THE FIRST 4 00000480 
CHARACTERS ARE REPRODUCED WITH THE FIRST PLOTTEDOO000490 
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POINT, AND THE LAST 4 WITH THE LAST POINT. SET 00000500 


TO BLANK ANY UNWANTED CHARACTERS, 00000510 

000 520 

1. ІТІТЕЕ AN ARRAY OF TWELVE 8-CHARACTER BCD WORDS. 00000530 
THE FIRST SIX WORDS WILL BE REPRODUCED AS THE 00000540 

FIRST LINE OF GRAPH TITLE» AND THE LAST SIX 00000550 


WORDS WILL FORM THE SECOND LINE. THE TITLE MUST 00000560 
INCLUDE THE USERS JOB IDENTIFICATIONe DIMENSION 00000570 


12 IN CALLING PROGRAM, AND SET TO BLANK ALL 00000580 
UNWANTED CHARACTERS. 00000590 
000 600 


8. EXSCACE X=SCALE (UNITS/INCH) AS POSITIVE FLOATING POINT 00000610 
VARIABLE WITH ONE FIGURE SIGNIFICANCEe SET TO 00000620 


ZERO FOR AUTO-SCALEo 00000630 

000 640 

9» YSCALE Y-SCALE (UNITS/INCH) AS POSITIVE FLOATING POINT 00000650 
VARIABLE WITH ONE FIGURE SIGNIFICANCEe SET TO 00000660 

ZERO FOR AUTO-SCALEo 00000670 

| 000 680 

10. ІХОР X-AXIS OFFSET FROM BOTTOM OF GRAPH IN INCHES, 00000690 
THIS MUST NOT EXCEED IHIGHs AND MUST NOT BE 00000700 

NEGATIVE. 00000710 

000 720 


11. IYRIGHT Y-AXIS OFFSET FROM LEFT OF GRAPH IN INCHES. THISOO000730 
MUST NOT EXCEED IWIDE, AND MUST NOT SE NEGATIVE 200000740 


000 750 
12. MODEXAX MODE OF X-AXIS OFFSET. SEE CODES BELOW, 00000760 
000 770 


13. MODEYAX MODE OF Y-AXIS OFFSET. THE CODES ARE AS FOLLOWS 00000780 
0 COMPUTED OFFSET, HOLDING ORIGIN ON 00000790 

GRAPH. THE CORRESPONDING IXUP OR 00000800 

IYRIGHT IS IGNORED 00000810 

1 COMPUTED OFFSET» WITH ORIGIN OFF THE 00000820 

GRAPH IF APPROPRIATE. THE CORRESPOND-00000830 

ING IXUP OR IYRIGHT IS IGNORED, USE | 00000840 


ONLY WITH AUTO-SCALE 00000850 

2 AXIS OFFSET AS SPECIFIED BY IXUP OR 00000860 

I YRIGHT, 00000870 

000 880 

14. IWIDE WIDTH OF GRAPH IN INCHES. THIS MUST NOT EXCEED 00000890 
NINE. ZERO WILL BE READ AS EIGHT INCHES, 00000900 

000 910 

15. IHIGH HEIGHT OF GRAPH IN INCHES. THIS MUST NOT EXCEED 00000920 
FIFTEEN. ZERO WILL BE READ AS EIGHT INCHES, 00000930 

000 940 

16. IGRID IF SET TO 1, A ONE INCH BY ONE INCH GRID WILL 00000950 
BE SUPERIMPOSED ON THE GRAPH. 00000960 

000 970 

17. LAST INDICATES TO CALLING PROGRAM WHETHER LAST PLOT 00000980 
WAS COMPLETED SUCCESSFULLY. THE CODES ARE 00000990 

0 LAST PLOT COMPLETED SUCCESSFULLY 00001000 

1 LAST PLOT NOT SUCCESSFUL 00001010 

2 LAST PLOT NOT SUCCESSFUL AND NO 00001020 
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ann ж 


9070 


9071 


1001 
1002 


1003 
1004 


1005 
1006 
1007 


FURTHER GRAPH OUTPUT WILL BE ATTEMPT-00001030 


ED UNTIL MOOCURV IS 
3 DRAW WAS ENTERED WITH MODCURV NOT 


EQUAL TO ONE OR ZERO WHILE THE ERROR 


LOCK-OUT WAS SET. 
THIS ARGUMENT MUST ALWAYS BE A NAME IN 
STATEMENT, NEVER A NUMBER, 


NOTE --— 
ALL ARGUMENTS FROM NUMBER 7 THRU, 
MODCURV IS EITHER 2 OR 3. HOWEVER» ARGUMENTS MUST NEVER BE 
OMITTED FROM THE CALLING STATEMENT, 
WHICH ARE THEN IRRELEVANT, 
VALUE IN THE CALL STATEMENT. NO VALUE IN THE CALLING PROGRAM 
WILL BE ALTERED BY THIS SUBROUTINE. 


NUMBER 16 ARE IGNORED WHEN 


REFERENCE --- 
THE BINARY TAPE FORMAT REQUIRED BY THE OFF-LINE PLOTTER IS 
DESCRIBED IN THE WRITEUP OF THE CDC 160A GRAPH PLOT PROGRAM 
(108812 8001) THE FORMAT REQUIRED BY THE COC 160 PROGRAM IS 
SIMILAR EXCEPT THAT THE INTERPOLATION ARGUMENT MUST BE ZERO. 


DIMENSION X1(900)» Y(900)» 
LTITLE(14)» 
JJTITLE(12) 

IPOINT = ITYPE 

CON ( ICONTRL = 400008) 

CONLICURV3 = 37773777202020208, ICURV4 = 0104000000000000B). 

REPLACE WITH DATA STATEMENT IN FORTRAN 62-30 


ITITLE(12)» JXTIT(12)» 
КАХІ5(5), ICURV(460) » 


JYTIT(12)» 
26КІ0(25)» ІСОМТ(1)» 


РОТ ІТЕЅТ = О IN DATA STATEMENT. 
IF (JTEST - 73546912) 9070»9071»9070 
ITEST = O 
JTEST = 73546912 
CONTINUE 


REMOVE ABOVE NONSENSE IN FORTRAN 63. 


CHECK PREVIOUS OPERATION OF ROUTINE, IF ANY, CODES ARE 
ITEST = 0 IF PREVIOUS GRAPH, IF ANY» COMPLETED 
ITEST = 1 IF PREVIOUS GRAPH NOT COMPLETED 
ITEST = 2 IF ERROR FOUND WHILE MODCURV WAS ONE» OR IF 


MODCURV WAS ILLEGAL, 

ІҒ(ІТЕ5Т - 2)1000»1001»91000 
IF (MODCURV ) 1003 100251003 
ITEST = O 
GO TO 1000 
IF (MODCURV - 1)10049 100221004 
LAST = 3 
RETURN 

SET UP ERROR RETURN ROUTINE. 
ІҒ(ІТЕ5Т)100991006»1009 
IF (MODCURV ) 10075100851007 
PRINT 1100 


ENTRY AT STATEMENT 1005, 


99 


NEXT ONE OR ZERO 


THE CALL 


IT IS MERELY THEIR VALUES 
ARGUMENTS MAY BE LISTED BY NAME OR 


00001040 
00001050 
00001060 
00001070 
00001080 
00001090 
00001100 
00001110 
00001120 
00001130 
00001140 
00001150 
00001160 
00001170 
00001180 
00001190 
00001200 
00001210 
00001220 
00001230 
00001240 
00001250 
00001260 
00001270 
00001280 
00001290 
00001300 
00001310 
00001320 
00001330 
00001340 
00001350 
00001360 
00001370 
00001380 
00001390 
00001400 
00001410 
00001420 
00001430 
00001440 
00001450 
00001460 
00001470 
00001480 
00001490 
00001500 
00001510 
00001520 
00001530 
00001540 
00001550 


1100 


1008 
1101 


1009 
1010 
1011 


1000 


100 


9000 
9100 


9001 
9002 
9003 
9101 


9004 
9005 
9102 


1020 


1021 
1102 


1022 
1023 
1024 
1103 


1 ONEe 


FORMAT ( 59H NO FURTHER GRAPH OUTPUT UNTIL MODCURV NEXT IS ZERO OR 

9/) 
ITEST = 2 
LAST = 2 
RETURN 
PRINT 1101 
FORMAT ( 30H THIS PLOT WILL NOT BE OUTPUT. 
LAST = 1 
RETURN 
IF(MODCURV - 2)010105,100851010 
IF(MODCURV - 3)10075101151007 
ITEST 3 O 
GO TO 1008 

CHECK LEGALITY OF INPUT ARGUMENTS, 

IF(NUMPTS - 2)1»2»2 
PRINT 100 
FORMAT (/» 
GO TO 1005 
IF (IPOINT)900029004,9001 
PRINT 9100 
FORMAT (/» 
GO TO 1005 
IF( IPOINT - 5)900259002,9000 
IF(NUMPTS - 30)3»3»9003 
PRINT 9101 
FORMAT (/» 
GO TO 1005 
IF(NUMPTS - 900)3,3,9005 
PRINT 9102 
FORMAT 1/» 
GO TO 1005 
IX = 1HX 
IY 1ΗΥ 
АМАХХ = -0.2E+100 
AMAXY = -0.2E+100 
AMINX = +0.2E+100 
AMINY = +0.2E+100 
DO 1020 I= 1l»sNUMPTS 
AMAXX = MAXIF(X( I) sAMAXX) 
AMAXY = MAXIF(Y( 1) sAMAXY) 
AMINX = MINIF(X(1) sAMINX) 
AMINY = MINIF(Y(1) sAMINY) | 
AMAXA = MAXIF(ABSF(AMAXX)» ABSF(AMAXY)» ABSFIAMINX)» ABSF(AMINY)) 
ІҒ(АМАХА - 1.0Е499)1022»1022»1021 
PRINT 1102 
FORMAT (/» 58H NO X OR Y VALUE MAY EXCEED 1.0E+99 IN ABSOLUTE MAGN 


7) 


32H NUMPTS MUST NOT BE LESS THAN 2. ) 


15H ILLEGAL ITYPEe ) 


46H NUMPTS MUST NOT EXCEED 30 FOR POINT PLOTTING. ) 


28H NUMPTS MUST NOT EXCEED 900, ) 


lITUDE, ) 


GO TO 1005 

ІҒ(АВОҒ(АМАХХ - AMINX) 
IF(ABSF(AMAXY - AMINY) 
PRINT 1103 
FORMAT (/» 
GO TO 1005 


- 1.0E-97)1023,1025»1025 
- 1.0E-97)1024,1025,1025 


38H ALL POINTS HAVE THE SAME COORDINATES, ) 
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00001560 
00001570 
00001580 
00001590 
00001600 
00001610 
00001620 
00001630 
00001640 
00001650 
00001660 
00001670 
00001680 
00001690 
00001700 
00001710 
00001720 
00001730 
00001740 
00001750 
00001760 
00001770 
00001780 
00001790 
00001800 
00001810 
00001820 
00001830 
00001840 
00001850 
00001860 
00001870 
00001880 
00001890 
00001900 
00001910 
00001920 
00001930 
00001940 
00001950 
00001960 
00001970 
00001980 
00001990 
00002000 
00002010 
00002020 
00002030 
00002040 
00002050 
00002060 
00002070 
00002080 


102 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 


104 


21 
22 
23 


24 
27 
28 


29 
30 
31 


32 


IFCITEST) 49794 


IF(MODCURV - 2)5,240,5,5 
IF(MODCURV - 3)6, 240,6 
PRINT 101 

FORMAT (/» 17H ILLEGAL MODCURV, 
GO TO 1005 

IF (MODCURV ) 62998 

IF (MODCURV - 1)69996 

IF (CIWIDE) 10011912 

ІТІТ = 5HIWIOE 

PRINT 102» ITIT» ITIT 
FORMAT (/» 9H ILLEGAL 5A55,29H, 
1 SH = Be »/) 
JWIOE = 8 

GO TO 14 

ІҒ(СІМІОЕ - 9)13,13,10 
JWIDE = IWIDE 

IF (I1HIGH) 15216917 

ITIT = SHIHIGH 

PRINT 102» ITIT» ITIT 
JHIGH = 8 

GO TO 19 

ІҒ(ІНІСН - 15)18,18,15 
JHIGH = IHIGH 

NODEXAX = MODEXAX 

IF (MODE XAX) 20927921 
ITIT= BHMODEXAX o 

PRINT 104, ITIT» IX 
FORMAT (/» 9H ILLEGAL 
1 Als ТНАХ = 0, »/) 
NODEXAX = O 

GO TO 27 

IF(MODEXAX - 2)275»22520 
IFCIXUP - JHIGH) 24924923 
ITIT = 8HIXUP, 

PRINT 1049 ITITs ІХ 
NODEXAX = O 

GO TO 27 

IF CIXUP) 23926926 

JXUP = IXUP 

NODEYAX = MODEYAX 

IF (MODEYAX ) 28 935929. 

ITI T*8HMODEYAXe 

PRINT 1049 ITITs IY 
NODEYAX = O 

GO TO 35 

ІҒ(МОРЕҮАХ – 2) 35,30, 28 
IF(IYRIGHT ~ ЈМІОЕ } 32,32,31 
ITIT = 8HIYRIGHT. 

PRINT 104, ITIT, IY 
NODEYAX = 0 

60 ТО 35 

IF(IYRIGHT) 31934934 


) 


GRAPH WILL BE PLOTTED WITH »А5» 


948» 32H GRAPH WILL BE PLOTTED WITH MODE, 
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00002090 
00002100 
00002110 
00002120 
00002130 
00002140 
00002150 
00002160 
00002170 
00002180 
00002190 
00002200 
00002210 
00002220 
00002230 
00002 240 
0000 2250 
00002260 
00002270 
00002280 
00002290 
00002300 
00002310 
00002220 
00002330 
00002340 
00002350 
00002360 
00002370 
00002380 
00002390 
00002400 
00002410 
00002420 
00002430 
00002440 
00002450 
00002460 
00002470 
00002480 
00002490 
00002500 
00002510 
00002520 
00002530 
00002540 
00002550 
00002560 
00002570 
00002580 
00002590 
00002600 
00002610 


C 
C 


34 JYRIGHT * IYRIGHT 00002620 
INITIALIZE PRIOR TO SCALING AND AXIS LOCATING, 00002630 
IFLAG = O FOR PASS WITH XDATAe IFLAG = 1 FOR PASS WITH YDATA, 00002640 
35 DO 2235 IOTA=1912 00002650 
2235 JJTITLECIOTA) - ITITLE(CIOTA) 00002660 
IFLAG = O 00002670 
BETA = Û0. 00002680 
SCALE = EXSCALE 00002690 
IAXIS = JYRIGHT 00002700 
MODE = NODEYAX 00002710 
ISIZE = JWIDE 00002720 
IXY = IX 00002730 
174 = IY 00002 740 
AMAX = AMAXX 00002750 
AMIN = AMINX 00002760 
GO TO 52 00002770 
50 IFLAG = 1 00002780 
BETA = O, 00002790 
SCALE = YSCALE 00002800 
IAXIS = JXUP 00002810 
MODE = NODEXAX 00002820 
ISIZE = JHIGH 00002830 
AMAX = AMAXY 00002840 
AMIN = AMINY 00002 850 
IXY = IY 00002860 
IYX = IX 00002870 
CHECK SCALE AND GO TO FIXED OR AUTO SCALE ROUTINES. 00002880 
52 IF(SCALE) 53,59,56 00002890 
53 PRINT 1145» IXY, IXY 00002900 
114 FORMAT (/» 9H ILLEGAL 5,A1539HSCALE4. GRAPH WILL BE PLOTTED wITH AU00002910 
1TO sAls TH-SCALEe »/) 00002920 
GO TO 59 00002930 
EXPRESS FIXED SCALE IN E FORMAT WITH ONE FIGURE SIGNIFICANCE. 00002940 
56 ІҒ(5СА(Е - 14,0E*499)57,53,53 00002950 
57 IF(SCALE - 140E-99)53,53,58 00002960 
58 CALL SCALEIT(SCALESISCAL10,FACTOR»1) 00002970 
SCALE = FACTOR*10,**ISCAL10 | 00002980 
CHECK AND COMPUTE AXIS LOCATION IF NECESSARYe FIXED SCALE 00002990 

CASE. ITAG = O IF ORIGIN ON GRAPH OR 1 IF IT IS SUPPRESSED. 00003000 . 
IF(MODE - 1)10325,103151030 000053010 
1030 ITAG =0 00003020 

GO TO 203 | 00003030 : 
1031 PRINT 1104 s IYXs IXYs» IXY 00003040 


1104 FORMAT (/s 
GO TO 59 
IF (ABSF(AMAX — AMIN) 
IF (ABSF (AMAX) 
IF (AMAX )10362103491037 
IAXIS = ISIZE/2 

GO TO 1030 

IAXIS = ISIZE 

GO TO 1030 


1032 
1033 
1039 
1034 


1036 
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5H MODE sA1ls24HAX MUST NOT BE 1 UNLESS 
1 (AUTO-SCALE)e¢ GRAPH WILL BE PLOTTED WITH AUTO »A1» 


TH=SCALE « 


- 1.0E-97)1033,1038,1038 
- 1.0E-97)1034,1039,1039 


sAls57THSCALE 15 000003050 ` 
9/) 00003060, 
00003070 
00003080 ՛ 


00003090 
00003100 
00003110 
00003120 
00003130 


00003140 ` 


i 





1037 


1038 
1040 


59 
60 


64 


65 
116 


68 


78 
80 
81 
82 


1982 
1983 


1984 
83 


84 
85 
86 


87 


1 


1 
2 


IAXIS = O 00003150 
GO TO 1030 00003160 
IF (SIGNF Cle sAMAX) ~ SIEGNF(1¢ sAMIN) )10409103951040 00003170 
ASIZE = ISIZE 00003180 
IAXIS - -АМІМ/(АМАХ - AMIN)*ASTZE +05 00003190 
GO TO 1030 00003200 
AUTO SCALE ROUTINE. 00003210 

IF (MODE - 1)60,64,69 00003220 
АМАХ - МАХІ1ІҒ(0»» АМАХ) 00003230 
АМІМ = МІМІР (О, » AMIN) 00003240 
IFCABSF(LAMAX — AMIN) - 1$40E-97)655,68,58 00003250 
PRINT 116, IXY, IXY, IYX 00003 260 
FORMAT (/» SH ALL »Al»47H VALUES EQUAL. AUTO SCALE POSSIBLE ONLY 00003270 
IF THE »A1»29H VALUES ARE NON-ZERO АМО МООЕ»А1» ТНАХ - 2» ) 00003280 
GO TO 1005 00003290 
ASIZE = ISIZE 00003300 
SCALE = (АМАХ - AMIN)/ASIZE 00003310 
GO TO 83 00003320 
IF (ABSF CAMAX — AMIN) — 1¢0E-97) 70974974 00003330 
IF CABSF CAMAX) - 1$.0E-97)71»75^, 74 00003340 
PRINT 1189 IXY 00003350 
FORMAT (/» 5H ALL »A1»38H VALUES ZERO. AUTO SCALE NOT POSSIBLE. )00003360 
GO TO 1005 7 00003370 
ІҒ(АМАХ - 1.0Е-97) 76975975 00003380 
ІҒ(І5ІГЕ - ІАХ1І5)77»76»977 00003290 
SCALE1 = Oo 00003400 
GO TO 78 00003410 
AXIS = IAXIS 00003420 
ASIZE = ISIZE 00003430 
SCALE] = AMAX/(ASIZE - AXIS) 00003440 
IF(AMIN + 1¢0E-97) 79979280 00003450 
ІҒ(ІАХІ5)81»80»81 00003%60 
SCALE2 = Os 00003470 
GO TO 82 00003480 
AXIS = TAXIS 00003490 
SCALE2 = -AMIN/AXIS 00003500 
IF(SCALE1l + SCALE2)1984,1982+1984 00003510 
PRINT 1983, IYX» IYX 00003520 
FORMAT (/» 56H NONE OF THE PLOT LIES ON THE GRAPH WITH THIS SPECIFO0003530 
IED »A1»47H-AXIS LOCATION. GRAPH WILL BE PLOTTED WITH ΜΟΡΕ.ΑΙ, 00003540 
THAX = Os 9/) 00003550 
МООЕ - 0 00003560 
GO TO 60 00003570 
SCALE = MAXIF(SCALEl» SCALE2) 00003580 
CALL SCALEIT(SCALE» ISCAL10,FACTOR»3) 00003590 
IF(FACTOR - 5.05)85»85,84 00003600 
FACTOR = 1 00003610 
15СА(10 = ISCAL10 + 1 00003620 
GO TO 90 00003630 
IFCFACTOR - 2402287587596 00003640 
FACTOR = 5 00003650 
GO TO 90 00003660 
IF(FACTOR - 1401)89,89,88 00003670 
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88 


89 
90 


91 
92 
93 


9% 
95 


96 
120 


99 


97 


200 


201 


202 


203 
204 


205 


FACTOR = 2 
GO TO 90 
FACTOR = 1 
SCALE s FACTOR*10, ** J[SCAL 10 
COMPUTE AXIS LOCATION IF NECESSARY, AUTO SCALE CASE, 
IF(MODE - 1)92,91,93 
ΙΕ(5ΙΟΝΕ( 1 νΑΜΑΧ) - 5ΙΟΝΕί 1 ο sAMIN) 192994992 
IAXIS = -AMIN/SCALE + 0,5 
ITAG = 0 
GO TO 203 
IF (AMAX 1959959200 
IAXIS = ISIZE 
BETA = -AMAX/SCALE 
ІҒ (ВЕТА – 1,Е+12) 99,99,96 
PRINT 120, IXY 


00003680 
00003690 
00003700 
00003710 
00003720 
00003730 
00003740 
00003750 
00003760 
00003770 
00003780 
00003790 
00003800 
00003810 
00003820 


FORMAT (/s 15H THE ORIGIN OF »Als 43H CANNOT BE OFFSET MORE THAN 100003830 


1.0€+12 INCHES. ) 


GO TO 1005 
IBETA = BETA 4 0.5 
BETA = -IBETA 


BETA IS THE NUMBER OF INCHES OF ORIGIN SUPPRESSION POSITIVE IF 


TRUE ORIGIN IS BELOW OR TO LEFT OF THE GRAPH. 
ΙΕ(ΒΕΤΑ { 1.) 9799799232 


ITAG = 1 
GO TO 203 
TAXIS = O 


BETA = AMIN/SCALE 

IF(BETA - 1leE+12)2019201996 
IBETA = BETA + 0.5 

BETA = IBETA 

ІЕ(ВЕТА - 1.)93,202,202 
ITAG = 1 


RELEASE RESULTS TO REMAINING PART OF PROGRAM. START SECOND 


PASS FOR Y VALUES IF NOT YET COMPUTED. 
IFC IFLAG) 20592045205 
SCALEX = SCALE 
IXFACT = FACTOR 


ІХ5С10 ғ І5СА110 
IXAXIS = TAXIS 
ITAGX = ITAG 
ISIZEX = ISIZE 
BETAX = BETA 

GO TO 50 

BETAY = BETA 
SCALEY = SCALE 
IYFACT = FACTOR 
1Ү5С10 = ISCAL1O 
IYAXIS = IAXIS 


NOW WRITE RECORDS, 
ITAGY = ITAG 
1512ЕҮ = ISIZE 


THIS COMPLETES CALCULATION OF SCALE FACTORS ETCe NOW GENERATE 


TAPE RECORDS» FIRST» THE SCALE FACTOR TITLES, 
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00003840 
00003850 
00003860 
00003870 
00003880 
00003890 
00003900 
00003910 
00003920 
00003930 
00003940 
00003950 
00003960 
00003970 
00003980 
00003990 
00004000 
00004010 
00004020 
00004030 
00004040 
00004050 
00004060 
00004070 
00004080 
00004090 
00004100 
00004110 
00004120 
00004130 
00004140 
00004150 
00004160 
00004170 
00004180 
00004190 
00004200 


206 


9206 


207 
208 


209 
210 


211 
212 
213 


214 
215 


216 


9074 


9075 


9080 
9081 


9084 
9085 


9090 


ӘХТІТ(1) = 8H1 х= 
JXTIT(2) = 8HSCALE = 
JXTIT(3) = ICODE(SCALEX) 
JXTITt4) = 8H UNITS/I 
JXTIT(5) = 8HNCH. 
JYTIT(1) = 8H1 Ya 
JYTIT(2) = 8HSCALE = 
JYTIT(3) * ICODE(SCALEY) 
JYTIT(4) x 8H UNITS/I 


JYTIT(5) s 8HNCH, 
DO 9206 126911 
ҮХТІТ(І) = 8H 
JYTIT(I) = 8H 
IF(ITAGX)2119211,207 
IF (BETAX) 20892085209 
JXTIT(7) = 8H ADD - 
GO TO 210 
JXTIT(7) ғ 8H ADD + 
JXTIT(8) * ICODE (BETAX*SCALEX) 
JXTIT(9) = BH UNITS T 
JXTIT(10)= 8НО ALL X 
JXTIT(11)* 8HVALUES, 
IF( ITAGY)21652165,212 
18 (88147) 21392139214 
JYTIT(7) = 8H ADD - 
GO TO 215 
JYTIT(7) = 8H ADD + 
JYTIT(8) = ICODE (BETAY*SCALEY) 
JYTIT(9) = 8H UNITS T 
JYTIT(10)=8HO ALL Y 
JYTIT(11)* 8HVALUES, 
ICONT(1) = ICONTRL + 4 
INSERT TITLE SIZE (02B) AHEAD OF MAIN TITLE RECORD. 
CALL ISHIFT6 (ITITLEs LTITLE) 
TEST FOR ALL BLANK TITLES. 
ICHECK = 8H 
DO 9075 125196 
IFCITITLEC(T)-ICHECK) 907499075,9074 
IFCITITLE(CT) ) 9080»9075»9080 
CONTINUE 
11711 = 1 
ICONT(1) = ICONT(1) - 1 
GO TO 9081 
111 = 0 
DO 9085 1=7,12 
IF (ITITLECI) - ICHECK) 90845,9085,9084^ 
IF (ITITLECI))9090,9085,9090 
CONT INUE 
IT2 = 1 
ICONT(1) = ICONT(1) - 1 
GO TO 9091 
IT2 = O 
NOW GENERATE AXES RECORDS, 
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00004210 
00004220 
00004230 
00004240 
00004250 
00004260 
00004270 
00004280 
00004290 
00004300 
00004310 
00004320 
00004330 
00004340 
00004350 
00004360 
00004370 
00004380 
00004390 
00004400 
00004410 
00004420 
00004430 
00004440 
00004450 
00004460 
00004470 
00004480 
00004490 
00004500 
00004510 
00004520 
00004530 
00004540 
00004550 
00004560 
00004570 
00004580 
00004590 
00004600 
00004610 
00004620 
00004630 
00004640 
00004650 
00004660 
00004670 
00004680 
00004690 
00004700 
00004710 
00004720 
00004730 


9091 


240 
9007 
9700 


9701 
242 


9261 
9242 


LFTMGN = 0*100 

IBOTMGN > 1%100 

ІН = LFTMGN 

JH IBOTMGN +IYAXIS*100 

ІН * ISIZEX*100 

IHL = LFTMGN * ISIZEX*100 - 107 
KAXIS(1) * IPACK12((IHs JH» LH» IHL) 
JHL = JH - 13 


IHL2 = -100 
IVH = (ISIZEX -IXAXIS - 1)*IXFACT 
IVH2 = -IXFACT 


KAXIS(2) = IPACK12 (JHL»o THL291VH»s IVH2) 

NH = ISIZEX 

ISH = 8H 14 

IV = LFTMGN * IXAXIS*100 

JV = IBOTMGN 

KAXIS(3) 7 IPACK12(NHs ISH9 IV9$ JV) 

LV = ISIZEY*100 

IVL = IV = 3 

JVL = IBOTMGN * ISIZEY*100 - 107 

JVL2 = -100 

KAXIS(4) = IPACK12(LV»IVL>»JVL,JVL2) 

IVV = (ISIZEY - IYAXIS - 1)*IYFACT 

IVV2 = -IYFACT 

INV = ISIZEY 

ISV = 8H 11 

KAXIS(5) = IPACK12(1VV» IVV2, INV»ISV) 
NOW GENERATE CURVES, 

SCX = 100,/SCALEX 

SCY = 100./SCALEY 

EXAXIS= IXAXIS*100 

YAXIS z IYAXIS*100 

ALFTMGN = LFTMGN 

BOTMGN = IBOTMGN 

SHIFTX =EXAXIS - BETAX*100. + ALFTMGN 

SHIFTY = YAXIS - BETAY*100. + BOTMGN 

EXSIZE= ISIZEX*100 + LFTMGN + 60 

SIZEX - LFTMGN - 60 

YSIZE = ISIZEY*100 + IBOTMGN + 70 

SIZEY = IBOTMGN - 70 

ICURV(1) = O 

IF (IPOINT 190109900 739010 

IF (XMODF (NUMPTSs2) 1970099701 99700 

ISWITCH = 1 

GO TO 242 

ISWITCH = 2 

INUM - (МОМРТ5 % 13/2 

DO 244 I=1»,INUM 

C] » X(2*I-1)*SCX + SHIFTX 

C2  Y(2*I-1)*5CY + SHIFTY 

IF CI-INUM) 241 99241 924] 

GO TO (92429241)sISWITCH 

mo s 
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00004740 
00004750 
00004760 
00004770 
00004780 
00004790 
00004800 
00004810 
00004820 
00004830 
00004840 
00004850 
00004860 
00004870 
00004880 
00004890 
00004900 
00004910 
00004920 
00004930 
00004940 
00004950 
00004960 
00004970 
00004980 
00004990 
00005000 
00005010 
00005020 
00005030 
00005040 
00005050 
00005060 
00005070 
00005080 
00005090 
00005100 
00005110 
00005120 
00005130 
00005140 
00005150 
00005160 
00005170 
00005180 
00005190 
00005200 
00005210 
00005220 
00005230 
00005240 
00005250 
00005260 


241 


9243 


244 
246 
9010 
247 
1260 
260 
261 


265 
9268 


9269 
9270 


9271 


9015 
270 


9020 
272 


оид 


1273 


1274 
1275 


C4 = C2 
GO TO 9243 
СЗ = Х(2%1)%5СХ + SHIFTX 
C4 к Ү(2%1)%5СҮ + SHIFTY 
Cl = MINIF(C1sEXSIZE) 
ІСІс МАХІҒ(СІ» 512ЕХ) 
C2 = МІМІЕ(С2» Ү512Е) 
ІС2: МАХІҒ(С2» 512ЕҮ) 
СЗ = MIN1F(C3»EXSIZE) 
IC3= MAXIF(C3, SIZEX) 
C4 = MINIF(C4s YSIZE) 
1Са= МАХІҒ( Са» SIZEY) 
ICURV(I*1) = IPACK12(1C1,1C291C39I1C4) 
II = INUM + 3 
CALL IPACKLI(LABEL» LABELI» IDUMMY) 
ICURV(II-1) = LABEL] 
ICURV(1I) = ICURV4 
IF(MODCURV - 1)247524759015 
CALL IREADY (IDUMMY) 
IF ( IDUMMY )5000 9126035000 
CALL IWRITE (ICONT» IDUMMY, 1) 
IF ( IOUMMY)500052605 5000 
CALL IWRITE (JXTIT9» IDUMMY 11) 
IF ( TIDUMMY )50009261,5000 
CALL IWRITE (JYTITs IDOUMMY 911) 
IF (IDUMMY ) 5000,265 5000 
IF(IT1)9269,9268,9269 
CALL IWRITE (LTITLE» IDUMMY, 7) 
IF ( IDUMMY ) 5000,9269, 5000 
IF (IT2)92719927099271 
CALL IWRITE (LTITLE(7)» IDUMMY, 7) 
IF t IDUMMY)50005,927155000 
CALL IWRITE (KAXIS» IDUMMY, 5) 
IF o IDUMMY )5000, 9015, 5000 
IF C IPOINT)902052705,9020 
CALL IWRITE (ICURV» IDUMMY, II) 
IF (o IDUMMY ) 5000, 9020, 5000 
IF (MODCURV - 1)27292729025 
IFCIGRID - 11)9025,273$,9025 
GENERATE GRID IF CALLED FOR. 
IX100 = ISIZEX*100 
IY100 s ISIZEY*100 
NEXT] z IBOTMGN 
NEXT2 = LFTMGN + IX100 
JGRID(1) = 0 
DO 1274 J=1,1992 
JGRID(J+1) с ІРАСК12 ((ІҒТМОМ» МЕХТІ» МЕХТ2» МЕХТ1) 
IF(NEXT1 - IBOTMGN - 1Ү100)1272»1275»1275 
NEXT1 = NEXT1 + 100 
JGRID(J*2) = IPACK12 (NEXT2» NEXT1» .LFTMGN» NEXT1) 
ІҒ(МЕХТІ - ІВОТМОМ - 1Ү100)127491276»91276 
ΝΕΧΤΙ = NEXT1 + 100 
JGRIO(J*2) - IPACK12 (NEXT2» NEXT1» NEXT2» NEXT1) 
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00005270 
00005280 
00005290 
00005300 
00005310 
00005320 
00005330 
00005340 
00005350 
00005 360 
00005370 
00005380 
00005390 
00005400 
00005410 
00005420 
00005430 
00005440 
00005450 
00005460 
00005470 
00005480 
00005490 
00005500 
00005510 
00005520 
00005530 
00005540 
00005550 
00005 560 
00005570 
00005580 
00005590 
00005600 
00005610 
00005620 
00005630 
00005640 
00005650 
00005660 
00005670 
00005680 
00005690 
00005700 
00005710 
00005720 
00005 730 
00005740 
00005750 
00005760 
00005770 
00005780 
00005790 


1276 


1277 


1278 


1279 
1280 
1281 


9025 


9030 


9031 
9032 
9033 
9034 


9035 


9036 


9037 


9038 


9039 


JGRID(J*3) = ICURV3 
JGRID(J+4) = ICURV4 

CALL IWRITE(JGRID»s IDUMMY » 
IF ( TIDUMMY )500091277,5000 
NEXT] = LFTMGN 

NEXT2 = IBOTMGN + IY100 
ОО 1279 Ј=1,11,2 


JGRID(J+1) = ТРАСК12 (NEXT1» 
ΙΕ(ΝΕΧΤΙ - LFTMGN - IX100)1278,1280,1280 


NEXT] = NEXT] + 100 


JGRID(J*2) = IPACK12 (МЕХТ1» 
IF(NEXT1 - LFTMGN - IX100)11279,1281,1281 


МЕХТІ = NEXT1 + 100 


JGRID(J*2) - ІРАСК12 (МЕХТ1» МЕХТ2» 


JGRID( J+3) 
JGRID(J+4) 


ICURV 3 
ICURV 4 


иы! 


CALL IWRITE (JGRID» IDUMMY, 


IF ( TIDUMMY ) 5000,9025,5000 
IF(IPOINT)9030s276,9030 


( 4+ ل 


J+4) 


ІВОТМСМ» 


МЕХТ2» 


NEXT1» 


NEXT1» 


NEXT1» 


GENERATE POINT PLOT RECORDS IF CALLED FOR, 


IOUT = O 


CALL IPACKL1 (LABEL» LABEL1» 


00 9050 Ігі1»МОМРТ5 
C] » X(I)*SCX + SHIFTX 
C2 - Ү(І)%5СҮ + SHIFTY 


ІЕ(СІ - ЕХ512Е)9031»9031»9034 
IF(C2 - YSIZE)9032,9032,9034 
ІЕ(СІ - 512ЕХ)9034»9033»9033 
IF(C2 - SIZEY)9034,9035,9035 


IOUT = IOUT +1 


GO TO 9050 
IC] = Cl 
IC2 = C2 


LABEL2) 


GO TO (9036 1903729038 »9039 59040) » IPOINT 


GENERATE CROSS, 
ICURV(2) * IPACK12 (IC1-5» 
ICURV(3) - ІРАСК12 (ІСІ » 
ICURV(^) = IPACK12 (IC1+5» 
GO TO 9041 

GENERATE PLUS. 
ІСОКУ(2) г ІРАСК12 (ІСІ » 
ICURV(3) г ІРАСК12 (ІСІ 9» 
ICURV(4) = IPACK12 (IC1+5» 
GO TO 9041 

GENERATE SQUARE, 
ICURV(2) = IPACK12 (1С1+4›, 
ICURV(3) = IPACK12 (IC1-4> 
ICURV(4) = IPACK12 (1C1+4» 
GO TO 9041 

GENERATE DIAMOND. 
ICURV(2) = IPACK12 (1C1+5» 
ICURV(3) ғ ІРАСК12 (ІСІ-5» 
ICURV(4) x IPACK12 (IC1*5» 


ІС2-5» 
IC2 9 
ІС2-5» 


ІС2-5, 
IC2 »9 
ІС2 9 


IC2-4, 
IC2445 
ІС2-4» 


Ico. 9 
IC2 > 
IC2 
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1С1+5, 
1C1-5» 
ІС1%5» 


ІСІ ,» 
IC1-5; 
1C1+5» 


ІСІЗ4» 
ІС1-4» 
1С1+4, 


ІСІ 9 
ІСІ » 
11459 


1C2+5) 
IC2 +5) 
IC2-5) 


1C2+5) 
ІС2 ) 
ІС2 ) 


1C2+4) 
IC2-4) 
IC2-4) 


1С2%5) 
1С2-5) 
ІС? ) 


МЕХТ2) 


IBOTMGN) 


NEXT2) 


00005800 
00005810 
00005820 
00005830 
00005840 
00005850 
00005860 
00005870 
00005880 
00005890 
00005900 
00005910 
00005920 
00005930 
00005940 
00005950 
00005960 
00005970 
00005980 
00005990 
00006000 
00006010 
00006020 
00006030 
00006040 
00006050 
00006060 
00006070 
00006080 
00006090 
06006100 
00006110 
00006120 
00006130 
00006140 
00006150 
00006160 
00006170 
00006180 
00006190 
00006200 
00006210 
00006220 
00006230 
00006240 
00006250 
00006260 
00006270 
00006280 
00006290 
00006300 
00006310 
00006320 


9040 
9041 
9042 


9043 
9044 


9045 
9046 
9050 


9048 
9104 


276 
277 
278 


130 


656 


279 


5000 
5001 


5002 
5100 


5670 


1МЕХ 


1 
IDUMMY = ITYP2( IDUMMY ) 


GO TO 9041 
GENERATE TRIANGLE. 
ICURV(2) » IPACK12 (IC1«*5» 
ICURV(3) = IPACK12 (1C1-5» 
ICURV(4) = ICURV(3) 
IF(I - NUMPTS)9043 99042 9043 
ICURV(5) = LABEL2 
GO TO 9046 
ІҒ(І - 1)9045»9044»9045 
ICURV(5) = LABEL1 
GO TO 9046 
ICURV(5) = ICURV3 
ICURV(6) = ICURV4 
CALL IWRITE (ICURVs IDUMMY, 6) 
IF  IDUMMY)5000590505, 5000 
CONTINUE 
IF(IOUT)9048,276590^8 
PRINT 9104, IOUT 
FORMAT (/» 1X» I2» 29H POINT(S) WERE OFF THE GRAPH, 
SET UP RETURN. 
IF (MODCURV) 27792785277 
IF (MODCURV - 3)279 92789279 
ITEST = О 
РКІМТ 1309 (.ҮӘТІТШЕ(І)»1:1912) 
FORMAT (/» 19H GRAPH TITLED o“. 


ІС2-2» ІСІ , 
IC2-3» IC1*5,» 


IC2 +6) 
1С2 3997 


96А8»/»919Х»;6А8, 
24H . oe HAS BEEN PLOTTED, »/»1H0) 


IF ( IDUMMY 156705656 55670 


LAST = 0 
RETURN 
ITEST = 1 
IDUMMY = ICLOCK( I DUMMY ) 
LAST = 0 
RETURN 
THESE ARE THE NORMAL RETURNS. 
NOW SET UP THE RETURN FOLLOWING A TAPE ERROR} 
IF (MODCURV ~ 1) 5001»5001»5002 


ІООММҮ = ITYPE1( IDUMMY ) 

GO TO 247 

PRINT 5100 

FORMAT (/» 36H TAPE ERROR IN WRITING GRAPH OUTPUT. 
IDUMMY = ITYPE1( IDUMMY ) 


GO TO 1007 
IDUMMY = ITYPE1( IDUMMY ) 
END 


SUBROUTINE IREADY (IDUMMY) 


SELECTS TAPE 8 (WILL LOOP UNTIL READY). WRITES EOF ОМ 8, 
MACHINE LANGUAGE WILL NOT BE NECESSARY IN FORTRAN 62-36 


EXF (520418) EXF7 (52000B), 
EXF7(00050B) SLJ (1RDY). 
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SELECT READ AND WAIT TAPE. 
EXIT ON CH 5 ACTIVE. 


00006330 
00006340 
00006350 
00006360 
00006370 
00006380 
00006390 
00006400 
00006410 
00006420 
00006430 
00006440 
00006450 
00006460 
00006470 
00006480 
00006490 
00006500 
00006510 
00006520 
00006530 
00006540 
00006550 
00006560 
00006570 
00006580 
00006590 
00006600 
00006610 
00006620 
00006630 
00006640 
00006650 
00006660 
00006670 
00006680 
00006690 
00006700 
00006710 
00006720 
00006730 
00006740 
00006750 
00006760 
00006770 
00006780 
00006790 
00006800 
00006810 
00006820 
00006830 
00006840 
00006850 


KAN 


1BUF 
1RDY 


2BUF 
1EOF 


1END 


1AGN 


1BUF 


2AGN 


1END 
2END 


LOC(ISIX = 
-EXF7(00061B). 


EXF7(52000B) 
LDA (IFIVE) 
ЕХҒЕ5(М)» 

EXF (520418) 
ΕΝΑ (0). 

ΕΧΕ (020008) 
EXF7(52001B) 


-EXF7(52000B). 


EXF7(52007B) 
ENA (IDUMMY ) 
INA (1) 
EXF5(N) 

ENA (0) 

EXF (62041B) 
ΕΧΕ (620418) 
ΕΧΕ (620038) 
ΕΧΕΤ( 620078) 
STA (IDUMMY)o 


END 


SUBROUTINE IWRITECISTART >» 
WRITE RECORD OF ІМОКО5» 


SLJ 
SAU 


(1МЕХ)» 
(1BUF )e 


EXF7(52000B), 


EXF 
SLJ 


(520068). 
(1END). 
SLJ (1E0F)e 
SAU (2BUF). 
STA (IFIVED, 
EXF7(52000B), 
EXF7(00061B). 
EXF7(62000B). 
EXF7(62000B). 
EXF7(62000B), 
ЕМА (10), 


I DUMMY » 
STARTING WITH ISTART. 


IF RECORD CORRECTLY WRITTEN» 


MACHINE 
6). 


ЕХҒ (620418) 
ЕХҒ (62041В) 
ΕΝΩ (1118). 
ENA (І5ТАКТ) 
SAL (1BUF) 
STA(ISIX) 
ЕМА(0) 
EXF7(62007B) 
EXF7(62003B) 
EXF7(62004B) 
EXF (62006B) 
QRS (3) 

EXF (62003B) 
STA (IDUMMY). 


END 


EXF7(62000B), 
EXF7(62000B), 


ІМА(1)» 
ADD(IWORDS), 
ΕΧΡ6(Ν). 
EXF7(62000B), 
SLJ (1END). 
ՏԼ (2AGN), 
SLJ (2bEND), 
EXF7(62000B), 
QJP1(1AGN), 
ЕМА (10). 


FUNCTION ITYP2 (IDUMMY ) 


TYPE WORD GRAPH, 
WILL NEED REWRITING IN FORTRAN 62-3, 
= 4513123015050000B ). 


CON(LC = 57B, 
LOC(ITWO Հ 2)» 
-ЕХЕ?Т( 000618)» 


EXF (620418) 
EXF (62041B) 
EXF (62003B) 
ENA (0) 

EXF7(62007B) 


EXF7(62000B), 
EXF7(62000B), 
ΕΧΕ (010008). 
EXF7(62000B), 
ΕΝΑ (10). 
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EXIT ON TAPE READY. 
TERMINATE 

BUFFER. 
SELECT AND WAIT TAPE 8e 
CLEAR Ae 
STOP CLOCK AND BACKSPACE 8, 
EXIT IF NOT AT LOAD POINT. 
WAIT TAPE 8. 
EXIT IF NO EOF, 
MOVE 

FORWARD 

OVER RECORD, 

CLEAR As WAIT CH 6. 
SELECT WRITE AND WAIT TAPE. 
SELECT AND WAIT TAPE 8. 
WRITE EOF AND WAIT. 
EXIT ON NO END OF TAPE. 


IWORDS) 
PUT IDUMMY = 


OTHERWISE SET NON-ZERO. 
LANGUAGE WILL NOT BE NECESSARY IN FORTRAN 62-36 


WAIT CH 66 

SELECT WRITE» WAIT ТАРЕ,» 
SELECT AND WAIT TAPE 8. 
SET COUNTER, 

STARTING ADDRESS, 
TERMINAL ADDRESS, 

BUFFER OUT. 

CLEAR Ae WAIT TAPE 8. 
EXIT IF NO END OF TAPE. 
EXIT IF NO PARITY ERROR, 
EXIT IF LENGTH ERROR. 
BACKSPACE AND WAIT. 

TRY WRITE 3 TIMESe 

WRITE EOFs NON ZERO Ae 
STORE RESPONSE, 


WAIT CH 6, 
SELECT AND 
SELECT AND 


WAIT TAPE, 
WAIT TAPE 8. 
WRITE EOFe START CLOCK. 
CLEAR Ae WAIT TAPE 8e 
EXIT IF NO END OF TAPE. 


00006860 
00006870 
00006880 
00006890 
00006900 
00006910 
00006920 
00006930 
00006940 
00006950 
00006960 
00006970 
00006980 
00006990 
00007000 
00007010 
00007020 
00007030 
00007040 
00907050 
00007060 
00007070 
00007080 
00007090 
00007100 
00007110 
00007120 
00007130 
00007140 
00007150 
00007160 
00007170 
00007180 
00007190 
00007200 
00007210 
00007227 
00007230 
00007240 
00007250 
00007260 
00007270 
00007280 
00007290 
00007300 
00007310 
00007320 
00007330 
00007340 
00007350 
00007360 
00007370 
00007380 


BAAN C 


1TYP 


1TYP 


STA (ITYP2), STORE RESPONSE. 


-EXF7(00021B). WAIT CH 2. 
EXF7(11141B) SLI (1TYP). EXIT IF UPPER CASE. 
ЕХҒ (211008) ENA (LC+1)e TYPE 
STA (ITWO) EXF2(LC) LOWER CASEe 

-EXF7(00021B). WAIT CH 2. 

EXF (21100B) ENA (М1+1), ТҮРЕ 
STA (ITWO) ЕХҒ2(М1)» GRAPH 
END 

FUNCTION ITYPE1 (IDUMMY) 


REWIND TAPE 8» REQUEST NEW ТАРЕ» 
WILL NEED REWRITING IN FORTRAN 62-3, 
CON(LC = 57Bs M1 = 4515112030242004B, M2 = 11033022040620318, 
M3 = 0401301520043342B). 
RSV(MESS = 3), 
LOC(LITWO = 2). 
-EXF7(00061B). 
EXF (62041B) 
EXF (62041B) 
EXF (62003B) 
ΕΧΕ (620078) 
ΕΧΕΤ(111418) 


AND WAIT TILL READY. 


WAIT СН 6. 

SELECT AND WAIT TAPE. 
SELECT AND WAIT TAPE 8e 
WRITE EOF AND WAIT. 
REWIND WITH INTERLOCK» 
EXIT IF UPPER CASE. 


EXF7(62000B), 
EXF7(62000B), 
EXF7(62000B), 
EXF7(00021B), 
SLJ (1TYP)o 


EXF (211008) ΕΝΑ (LC+1)e TYPE 

STA (I TWO) EXF2(LC). LOWER CASE. 
LDA (M1) STA (MESS), 

LDA (M2) STA (МЕ5541), 

LDA (M3) STA (MESS*2), 
-EXF7(00021B). WAIT CH 2 
ЕХҒ (211008) ENA (MESS+3), TYPE 

STA (ITWO) EXF2(MESS)e MESSAGE, 
-ЕХҒТ(000618)» WAIT СН 6. 
-EXF7(62000B). WAIT TAPE. 


EXF (620418) 
EXF (01000B). 
END 


EXF7(62000B), SELECT AND WAIT TAPE 80 


START CLOCK. 


FUNCTION ICODE (ANUMBER) 
CODES ABSOLUTE VALUE OF A FLOATING POINT NUMBER (BETWEEN 
1.OE-100 AND 1.0E+100) INTO 8-CHARACTER BCD WORD OF THE FORM 
1e23E+45. ICODE = 8HO+00E+00 IF MAGNITUDE OUT OF RANGE, 
CHECK AVAILABILITY OF LIBRARY FUNCTIONS IN FORTRAN 62-3, 
DIMENSION II(8) 
BNUMBER = ABSF(ANUMBER) 
IF(BNUMBER - 14,0E*100)75656 
IF(BNUMBER - 140E-100)65652 
ІСООЕ : 8Н06.00Е%00 
RETURN 
THIS IS ERROR EXIT. 
CALL SCALEIT (BNUMBER, ISCAL10; 
ISIGNSC » XSIGNF(19»1SCAL10) 
ISCAL10 * XABSF(ISCAL10) 
IFACT = FACTOR*100.001 


FACTOR, 3) 


ՅԼ 


00007390 
00007400 
00007410 
00007420 
00007430 
00007440 
00007450 
00007460 
00007470 
00007480 
00007490 
00007500 
00007510 
00007520 
00007530 
00007540 
00007550 
00007560 
00007570 
00007580 
00007590 


WAIT CH 200007600 


00007610 
00007620 
00007630 
00007640 
00007650 
00007660 
00007670 
00007680 
00007690 
00007700 
00007710 
00007720 
00007730 
00007740 
00007750 
00007760 
00007770 
00007780 
00007790 
00007800 
00007810 
00007820 
00007830 
00007840 
00007850 
00007860 
00007870 
00007880 
00007890 
00007900 
00007910 


ο + 


* rx rn СУ С\ б 


11(8) = XMODF(ISCAL10,10) 00007920 
II(7)*15CAL10/10 00007930 
IF (ISIGNSC) 49393 00007940 
11(6) = 8H + 00007950 
GO TO 5 00007960 
11(6) = 8H - 00007970 
11(5) = 8H E 00007980 
11(%) = XMODF(IFACT, 10) 00007990 
11(3) = (XMODF(IFACT»100))/10 00008000 
11(2) = 8H : 00008010 
11(1) = IFACT/100 00008020 
CALL IPACK (IIs IPACKED) 00008030 
ICODE = IPACKED 00008040 
RETURN 00008050 
END 00008060 
00008070 

SUBROUTINE SCALEIT (ANUMBER»ISCAL10+ FACTOR »MODE) 00008080 
FINDS FACTOR (BETWEEN 1+0 AND 9.99...) AND SCALE OF 10 AS 00008090 
DEFINED BY ANUMBER = FACTOR*10,'**ISCAL10, 00008100 
MODE IS THE NUMBER OF SIGNIFICANT FIGURES REQUIRED. THIS MUST 00008110 

BE BETWEEN 1 AND 10 OR IT WILL BE PUT EQUAL TO SIX, 00008120 
CHECK AVAILABILITY OF LOG10F IN FORTRAN 62-3. 00008130 
1SCAL10=LOG10F ( ANUMBER ) 00008140 
FACTOR = ANUMBER/10.**ISCAL10 00008150 
IF(FACTOR ~ 0,1)1,2,2 00008160 
FACTOR * FACTOR*100, 00008170 
ISCAL10 = ISCAL10 - 2 00008180 
GO TO 8 00008190 
IF(FACTOR – 1,0)3,8,4 00008200 
FACTOR - FACTOR*10, 00008210 
ISCAL10 = ISCAL10 - 1 00008220 
GO TO 8 00008230. 
IF (FACTOR - 100.0)69,595 00008240 
FACTOR = FACTOR/1000 00008250 
ISCAL10 = ISCAL10 + 2 00008260 
GO TO 8 00008270 
ІЕ(ҒАСТОВ - 10.0)8»7»7 00008280 
FACTOR = FACTOR/10. 00008290 
ISCAL10 = ISCAL10 + 1 00008300 
IF (MODE)9,9,10 00008310 
MODE = 6 00008320 
GO TO 11 | Е 00008330 
IF (MODE ~- 10)1151159 00008340 
IFACTOR = FACTOR£*10,**(MODE - 1) 4 0.5 00008350 
FACTOR = IFACTOR 00008360 
FACTOR = FACTOR/10.**(MODE - 1) 00008370 
IF(FACTOR – 10, )13,12,12 | 00008380 
FACTOR = le 00008390 
ISCAL10 = ISCAL10 + 1 00008400 
RETURN 00008410 
END 00008420 
00008430 

SUBROUTINE IPACK (IIs IPACKED) 00008440 


2) 


ROM 


1NEX 


2NEX 


* ^ 


1NEX 


2NEX 


TAKES 8 SIX-BIT WORDS AND PACKS THEM LEFT TO RIGHT 
IN IPACKED. IF WORD IS ZERO» 12B IS SUBSTITUTED. 
CONVERT TO CODAP FOR FORTRAN 62-36 

CON(IZERO к 128). 


SIU1(ISAVE) ЕМ11(8). 
LDA1(11) AJP1(2NEX)o 
LDA (IZERO). 

LRS (6) INI1(-2). 
15К1(-1) SLJ (1NEX), 
STQ (IPACKED) — LIUICISAVE) 
END 


SUBROUTINE ISHIFT6 (ITITLE» LTITLE) 
INSERTS 02B AHEAD OF 6-WORD TITLE RECORD, 
WILL HAVE TO BE CONVERTED TO СОРАР IN FORTRAN 62-3, 
WATCH ARRAY INDEXING IN FORTRAN 62-3. 
CON ( IBLANK = 2020202020202020B)0 


SIUICISAVE) ЕМІ1(1)» SAVE INDEX» SET COUNT, 
ΕΝΑ (2). ENTER O2Be 
LOQICITITLE) LLS (42). PERFORM 
STA1(LTITLE) LLS (6). SHIFTING 

15К1(6) SLI (1МЕХ)» CHECK IF COMPLETE. 
LDQ (IBLANK) LLS (42). COMPLETE LAST WORD. 
ENT1(7) STAL(LTITLE). STORE LAST. 

ENA (2), REPEAT 

LDG1(1TITLE) ԼԼՏ (42). FOR 

ІМІ1(1) STAL(LTITLE). SECOND 

INI1(-1) ԼԼՏ (6). TITLE 

ISK1(12) SLJ (2ΝΕΧ)ο LINE. 

LDQ (IBLANK) LLS (42). 

ENT1(14) STA1(LTITLE). 

LIUICISAVE), RESTORE INDEX, 
END 


FUNCTION IPACK12 (10МЕ»12913»14) 
PACKS FOUR 12-BIT WORDS INTO ONE 48-BIT WORD. 
WILL REQUIRE CONVERSION TO CODAP IN FORTRAN 62-3. 


LDA( IONE) LDQ(I2). 
QLS (36) LLS(12)0 
(00(13) QLS( 36). 
LLSt12) LDQ( 14). 
QLS(36) LLS(12)0 
STA(IPACK12). 

END 


SUBROUTINE IPACKL1 (LABEL» LABEL1» LABEL2) 
PACKS TWO 4-CHARACTER LABELS, 
USE DECODE/ENCODE IN FORTRAN 62-3. 
CON(IFLAG = 37773777B)e 
LDA (IFLAG) LDQ (LABEL). 
LLS(24) STA (LABEL1), 
LDA (IFLAG) LLS(24). 
STA (LABEL2)o 
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00008450 
00008460 
00008470 
00008480 
00008490 
00008500 
00008510 
00008520 
00008 530 
00008540 
00008550 
00008560 
00008570 
00008580 
00008590 
00008600 
00008610 
00008620 
00008630 
00008640 
00008650 
00008660 
00008670 
00008680 
00008690 
00008700 
00008710 
00008720 
00008730 
00008740 
00008750 
00008760 
00008770 
00008780 
00008790 
00008800 
00008810 
00008820 
00008830 
00008840 
00008850 
00008860 
00008870 
00008880 
00008890 
00008900 
00008910 
00008920 
00008930 
00008940 
00008950 
00008960 
00008970 


END 


FUNCTION ICLOCK(IOUMMY) 
ЕХҒ (010008). 

ЕМО 

END 


START CLOCK, 
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00008980 
00008990 
00009000 
00009010 
00009020 
00009030 


907 





APPENDIX III 


1. Derivations relating ՉԸ with да Жапа д with ad [4]. 
òX OX дү ay 
The equation С = 27 tank {24 ее) 
on rearrangement gives Հոմ = LM E ‚ 
C E 


The power series representation of the inverse hyperbolic tangent for 


values of (21)? « / is 
-//2πεὶ Fa Zn C)? шй 
tanh “> = (4 7 ON + £ (288 M ond 


so that to a good degree of accuracy it can be written 


x 2 T CE - «(2-9 (¥ 7 


and Ғы 
E pS LESS ее” 
Then dd 2 2 ç b 7 
2 «2 t UDS + RS AIG) 
Since the depth may be considered as being a function of C (а = rc], 


and C = G(X,Y), then by the chain rule 


dd | M ὃς 4 Dd ὃς 


х De 2X Τ.Θ 


also since d is an explicit function of x and y 


дЫ Է. 24 


ee յռ [ροι αἲ Ր NE 
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2. Derivations of the partial derivatives of the depth function with 


respect to X and Y. 
The derivatives may be computed directly from the equation for the 


quadric surface 


d = A, + A, X + Á; V + A XY + Z X 4 A 


so that 
dd . A. AY uiii 
Бу 7 E 7 5 
44 2, ἑΑ.Χ Zur 
and 
2 2 
— = 2 od — Ж 
оу уын 2-4. 
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APPENDIX IV 
Summary of refraction and shoaling relationships for 
intermediate depths [5] (6 |. 

Тһе following discussion applies to waves of small steepness where 
the deep-water wave height divided by the deep-water wave length is less 
than .005 (Ho/Lo < .005). In all cases the subscript zero refers to deep- 
water parameters. 

The wave velocity depends upon wave length and upon the depth of 


waters: 


. 25 4 (294 
С = 21 fanh (E 


where d is the depth, and T is the period of the wave. 

Waves of a certain period curve as they approach the shore from deep 
water until, theoretically, they are perpendicular to the beach in zero 
depth of water. For any change in depth, Snell's law determines the cur- 
vature of the ray. It must intersect a contour at an angle determined by 
Snell's law for the successive changes in depth. The tangent to the wave 
ray must make an angle, Ж, with a perpendicular to the contour at the 
point where the ray intersects the contour. The ray must curve with the 
change in depth so that Snell's law is satisfied at a discrete set of 
points given by the intersection of the ray with a set of contours. As 
shown in Figure IV, the wave ray crosses the contour corresponding to the 


wave speed С The tangent to the wave ray makes an angle X with a line 


1" 
drawn perpendicular to the smoothed contour. Since the wave ray is con- 
tinuously changing direction, it must make a new angle, X, with the 


perpendicular to the contour corresponding to the wave speed C when it 


2? 
reaches that contour. The change in angle is Δα. Then at the two con- 


tours corresponding to wave speeds, C, and SE Snell's law holds «Since 
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the wave crests intersect the contours at the correct angles. The impor- 
tant point is that the two ray tangents are connected by an arc of a 
circle which determines the exact path of the wave ray from point A to 
point B. The iteration procedure described in the text is used to arrive 
at this result. 

The assumption behind the wave height calculation is that for steady 
state conditions energy does not flow across orthogonals and that none is 
destroyed by friction. Therefore, the power between orthogonals is as- 
sumed to remain constant. The mean wave energy per unit surface area 


equals: 


£= g ea ° 


where g Ls the density of the water and H is the wave height. According 
to wave theory only a fraction of the wave energy is carried forward with 
the wave form at the speed C. Then the mean power transmitted between 


orthogonals equals: 


p=nEC di 


where n is the fraction of energy carried forward and dl is distance be- 
tween orthogonals. The numerical value of n approaches % in deep water 
and approaches one in shallow water. By equating the energy in the deep 
water to that in the shallow water, the ratio is formed: 


L£L n / 


E “2 


տամ. 


μέ / 
Դ C/Co 20/42, 


where the terms аге defined above. This can be written: 
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where K is the refraction coefficient and is equal to 41, /al. The 


term under the square root sign is termed the shoaling factor, H.: 
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